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On the arcminute angular scales probed by Planck, the CMB anisotropics are gently perturbed by gravitational lensing. Here we present a detailed 
study of this effect, detecting lensing independently in the 100, 143, and 217 GHz frequency bands with an overall significance of greater than 
25cr. We use the temperature-gradient correlations induced by lensing to reconstruct a (noisy) map of the CMB lensing potential, which provides 
an integrated measure of the mass distribution back to the CMB last- scattering surface. Our lensing potential map is significantly correlated with 
other tracers of mass, a fact which we demonstrate using several representative tracers of large-scale structure. We estimate the power spectrum 
• • of the lensing potential, finding generally good agreement with expectations from the best-fitting ACDM model for the Planck temperature power 
. ^ spectrum, showing that this measurement at z = 1100 correctly predicts the properties of the lower-redshift, later- time structures which source 
k>j the lensing potential. When combined with the temperature power spectrum, our measurement provides degeneracy-breaking power for parameter 
r\ constraints; it improves CMB -alone constraints on curvature by a factor of two and also partly breaks the degeneracy between the amplitude of the 
^ primordial perturbation power spectrum and the optical depth to reionization, allowing a measurement of the optical depth to reionization which 
is independent of large-scale polarization data. Discarding scale information, our measurement corresponds to a 4% constraint on the amplitude 
of the lensing potential power spectrum, or a 2% constraint on the RMS amplitude of matter fluctuations at z ~ 2. 



This paper, one of a set of papers associated with the 2013 release of data from the Planck" mission (Planck Collaboration I 
2013), describes our reconstruction of the CMB lensing potential based 15 months of data, estimation of the lensing potential power 
spectrum, and a first set of associated science results. 

^ Corresponding author: K. Benabed <benabed@iap . £r> 

* Corresponding author: D. Hanson <dhanson@physics .mcgill . ca> 

* Planck (http://www.esa.int/Planck) is a project of the European Space Agency (ESA) with instruments provided by two scientific 
consortia funded by ESA member states (in particular the lead countries Erance and Italy), with contributions from NASA (USA) and telescope 
reflectors provided by a collaboration between ESA and a scientific consortium led and funded by Denmark. 
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1. Introduction 

When Blanchard and Schneider first considered the eff'ect of 
gravitational lensing on the cosmic microwave background 
(CMB) anisotropics in 1987, they wrote with guarded optimism 
that although "such an observation is far from present possibil- 
ities [...] such an eff'ect will not be impossible to find and to 
identify in the future." (Blanchard & Schneider 1987). In the 
proceeding years, and with the emergence of the concordance 
ACDM cosmology, a standard theoretical picture has emerged, 
in which the large-scale, linear structures of the Universe which 
intercede between ourselves and the CMB last- scattering sur- 
face induce small but coherent (Cole & Efstathiou 1989) de- 
flections of the observed CMB temperature and polarisation 
anisotropics, with a typical magnitude of T . These deflec- 
tions blur the acoustic peaks (Seljak 1996), generate small-scale 
power (Linder 1990; Metcalf & Silk 1997), non-Gaussianity 
(Bernardeau 1997), and convert a portion of the dominant E- 
mode polarisation to 5-mode (Zaldarriaga & Seljak 1998). 
Gravitational lensing of the CMB is both a nuisance, in that it 
obscures the primordial fluctuations (Knox & Song 2002), as 
well as a potentially useful source of information; the charac- 
teristic signatures of lensing provide a measure of the distri- 
bution of mass in the Universe at intermediate redshifts (typi- 
cally 0.1<z<5). In the ACDM framework, there exist accurate 
methods to calculate the efl'ects of lensing on the CMB power 
spectra (Challinor & Lewis 2005), as well as optimal estimators 
for the distinct statistical signatures of lensing (Hu & Okamoto 
2002; Hirata & Seljak 2003a). 

In recent years there have been a number of increasingly sen- 
sitive experimental measurements of CMB lensing. Lensing has 
been measured in the data of the WMAP satellite both in cross- 
correlation with large-scale-structure probed by galaxy surveys 
(Hirata et al. 2004; Smith et al. 2007; Hirata et al. 2008; Feng 
et al. 2012a), as well as internally at lower signal-to-noise (Smidt 
et al. 2011; Feng et al. 2012b). The current generation of low- 
noise, high-resolution ground-based experiments has done even 
better; the Atacama Cosmology Telescope (ACT) has provided 
an internal detection of lensing at 4.6(T (Das et al. 2011, 2013), 
and the South Pole Telescope detects lensing at 8.1 cr in the tem- 
perature power spectrum, and 6.3(T from a direct reconstruction 
of the lensing potential (Keisler et al. 2011; van Engelen et al. 
2012; Story et al. 2012). Significant measurements of the cor- 
relation between the reconstructed lensing potential and other 
tracers of large-scale structure have also been observed (Bleem 
et al. 2012; Sherwin et al. 2012). 

Planck enters this field with unique full-sky, multi-frequency 
coverage. Nominal map noise levels for the first data release (ap- 
proximately 105, 45, and 60yuKarcmin for the three CMB chan- 
nels at 100, 143, and 217 GHz respectively) are approximately 
five times lower than those of WMAP (or twenty five times lower 
in power), and the Planck beams (approximately 10', 1' and 5' 
at 100, 143 and 217 GHz), are smaU enough to probe the 2.'4 
deflections typical of lensing. Full sky coverage is particularly 
beneficial for the statistical analysis of lensing eff'ects, as much 
of the "noise" in temperature lens reconstruction comes from 
CMB fluctuations themselves, which can only be beaten down 
by averaging over many modes. 

Lensing performs a remapping of the CMB fluctuations, 
such that the observed temperature anisotropy in direction h 
is given in terms of the unlensed, "primordial" temperature 



anisotropy as (e.g. Lewis & Challinor 2006) 

T{h) = T''''\h -h V0(w)), 

= T^'^'ih) + 2 ^'m^iT{h) + 0{^\ (1) 

i 

where 0(w) is the CMB lensing potential, defined by 

m = -2 r dx l^,^:~f^ ^(xn; m -X\ (2) 

Jo fK(X*)fK(X) 

Here;^ is conformal distance (with;^:, ^ 14000 Mpc) denoting 
the distance to the CMB last- scattering surface) and ^(xn,r]) 
is the gravitational potential at conformal distance x along the 
direction h at conformal time rj (the conformal time today is de- 
noted as tjq). The angular-diameter distance fxix) depends on 
the curvature of the Universe, and is given by 

( R-^'^ ^m{K^l^x) for > (closed), 
fKix)-\x for /r = (flat), (3) 

\\K\-^'^ sinh(|^|i/2^) for ^ < (open). 

The lensing potential is a measure of the integrated mass distri- 
bution back to the last-scattering surface. It contains information 
on both the gravitational potentials ^ To first order, its eff'ect on 
the CMB is to introduce a correlation between the lensed tem- 
perature and the gradient of the unlensed temperature, a property 
which can be exploited to make a (noisy) reconstruction of the 
lensing potential itself. 

In Fig. 1 we plot the noise power spectrum for recon- 
struction of the lensing potential using the three Planck frequen- 
cies which are most sensitive to the CMB anisotropics on the 
arcminute angular scales at which lensing eff'ects become ap- 
parent. The angular size of the Planck beams (5' FWHM and 
greater) does not allow a high signal-to-noise {S IN) reconstruc- 
tion of the lensing potential for any individual mode (our high- 
est S IN ratio on an individual mode is approximately 2/3 for the 
143 and 217 GHz channels, or 3/4 for a minimum- variance com- 
bination of both channels), however with full-sky coverage the 
large number of modes which are probed provides considerable 
statistical power. To provide a feeling for the statistical weight of 
diff'erent regions of the lensing measurement, in Fig. 2 we plot 
(forecasted) contributions to the total detection significance for 
the potential power spectrum C^j^ as a function of lensing mul- 
tipole L. In addition to the power spectrum of the lensing po- 
tential, there is tremendous statistical power in cross-correlation 
of the Planck lensing potential with other tracers of the matter 
distribution. In Fig. 2 we also plot forecasted S IN contributions 
for several representative tracers. 

This paper describes the production, characterization, and 
first science results for two Planck-dtrntd lensing products: 

(I) A map of the CMB lensing potential ^{n) over a large 
fraction of the sky (approximately 70%). This repre- 
sents an integrated measure of mass in the entire visible 
Universe, with a peak sensitivity to redshifts of z ~ 2. 
At the resolution of Planck, this map provides an esti- 
mate of the lensing potential down to angular scales of 
5' at L - 2048, corresponding to structures on the order 
of 3 Mpc in size at z = 2. 
(II) An estimate of the lensing potential power spec- 
trum C^j^ and an associated likelihood, which is 
used in the cosmological parameter analysis of 
Planck CoUaboration XVI (20 1 3). Our likelihood is 
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Fig. 1. Sky-averaged lens reconstruction noise levels for the 100, 
143, and 217 GHz Planck channels (red, green, and blue solid, 
respectively), as well as for experiments that are cosmic- variance 
limited to a maximum multipole ^max = 1000, 1500, and 1750 
(upper to lower solid grey lines). A fiducial ACDM lensing po- 
tential using best-fit parameters to the temperature power spec- 
trum from Planck Collaboration XVI (2013) is shown in dashed 
black. The noise level for a minimum-variance ("MV") combi- 
nation of 143 -h 217 GHz is shown in black (the gain from adding 
100 GHz is negligible). 
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Fig. 2. Overview of forecasted contributions to the detection sig- 
nificance as a function of lensing multipole L for the C^^ power 
spectrum (solid black), as well as for several other mass tracers, 
at the noise levels of our MV lens reconstruction. Our measure- 
ment of the power spectrum C^^ is presented in Sect. 6, The 
ISW-0 correlation believed to be induced by dark energy is stud- 
ied in Sect. 6.2. The NVSS-0 correlation is studied (along with 
other Galaxy correlations) in Sect. 6.3. The CIB-0 prediction 
(dashed cyan) uses the linear SSED model of Hall et al. (2010), 
assuming no noise or foreground contamination. A full analy- 
sis and interpretation of the CIB-0 correlation is performed in 
Planck Collaboration XVIII (2013). 



based on the lensing multipole range 40 < L < 400. 
This multipole range (highlighted as a dark grey band 
in Fig. 2), was chosen as the range in which Planck 
has the greatest sensitivity to lensing power, encap- 
sulating over 90% of the anticipated signal-to-noise, 
while conservatively avoiding the low-L multipoles 
where mean-field corrections due to survey anisotropy 
(discussed in Appendix C) are large, and the high-L 
multipoles where there are large corrections to the power 
spectra from Gaussian (disconnected) noise bias. Distilled 
to a single amplitude, our likelihood corresponds to a 
4% measurement of the amplitude of the fiducial ACDM 
lensing power spectrum, or a 2% measurement of the 
amplitude of the matter fluctuations (neglecting parameter 
degeneracies). 

Our efl'orts to validate these products are aided by the fre- 
quency coverage of the three Planck channels that we employ, 
which span a wide range of foreground, beam, and noise prop- 
erties. For the mask levels that we use, the root-mean- squared 
(RMS) foreground contamination predicted by the Planck sky 
model (Delabrouille et al. 2012) has an amplitude of 14, 22, 
and 70yuK at 100, 143, and 217 GHz, which can be compared 
to a CMB RMS for the Planck best-fitting ACDM power spec- 
trum of approximately llOyuK. The dominant foreground com- 
ponent at all three CMB frequencies is dust emission, both from 
our Galaxy as well as the cosmic infrared background (CIB), 
although at 100 GHz free-free emission is thought to consti- 
tute approximately 15% of the foreground RMS. Contamination 
from the thermal Sunyaev-Zeldovich (tSZ) eff'ect is a potential 
worry at 100 and 143 GHz, but negligible at 217 GHz (Sunyaev 
& Zeldovich 1980). On the instrumental side, these frequency 
channels also span a wide range of beam asymmetry, with typi- 
cal ellipticities of 19%, 4%, and 18% at 100, 143, and 217 GHz. 
The magnitude of correlated noise on small scales (due to de- 
convolution of the bolometer time response) also varies signifi- 
cantly. The ratio of the noise power (before beam deconvolution) 
at ^ = 1500 to that at ^ = 500 is a factor of 1.5, 1.1, and 1.0 at 
100, 143, and 217 GHz. The agreement of lens reconstructions 
based on combinations of these three channels allows a powerful 
suite of consistency tests for both foreground and instrumental 
biases. We will further validate the robustness of our result to 
foreground contamination using the component separated maps 
from the Planck consortium (Planck Collaboration XII 2013). 

At face value, the 4% measurement of C^j^ in our fiducial 
likelihood corresponds to a 25cr detection of gravitational lens- 
ing eff'ects. In fact, a significant fraction (approximately 25% of 
our error bar) is due to sample variance of the lenses themselves, 
and so the actual "detection" of lensing eff'ects (under the null 
hypothesis of no lensing) is significantly higher. We have also 
been conservative in terms of mask and multipole range in the 
construction of our fiducial lensing likelihood. As we will show 
in Sect. 7.1, we obtain consistent results on sky fractions larger 
than our fiducial 70% sky mask. 

The Planck lensing potential is part of a significant shift for 
CMB lensing science from the detection regime to that of preci- 
sion cosmological probe. The NVSS quasar catalogue, for exam- 
ple, has been a focus of previous lensing cross-correlation stud- 
ies with WMAP (Hirata et al. 2004; Smith et al. 2007; Hirata 
et al. 2008), where evidence for cross-correlation was found at 
approximately 3.5(T. As we will see in Sect. 6.3, the significance 
for this correlation with Planck is now 20(T. Notably, this is less 
than the significance with which lensing may be detected inter- 
nally with Planck. The lensing potential measured by Planck 



3 



Planck Collaboration: Gravitational lensing by large-scale structures with Planck 



now has sufficient signal-to-noise that shot noise of the NVSS 
quasar catalogue is the limiting source of noise in the cross- 
correlation. 

The majority of this paper is dedicated to the production and 
testing of the Planck lensing map and power spectrum estimate. 
Our focus here is on extracting the non-Gaussian signatures of 
lensing, although we note that lensing effects are also apparent 
at high significance (lOcr) as a smoothing eff'ect in the Planck 
temperature power spectra (Planck Collaboration XV 2013). We 
begin in Sect. 2, where we describe and motivate our method- 
ology for producing unbiased estimates of the lensing potential 
and its power spectrum. The Planck maps and data cuts that are 
used for this purpose are described in Sect. 3, and the simulations 
that we use to characterize our reconstruction and its uncertain- 
ties are described in Sect. 4. In Sect. 5 we give an overview of 
our error budget, and discuss the various sources of systematic 
and statistical uncertainty for our lensing estimates. In Sect. 6 
we present our main results: the first Planck lensing map and 
a corresponding estimate of the lensing potential power spec- 
trum. The likelihood based on this power spectrum is combined 
with the likelihood for the temperature anisotropy power spec- 
trum (Planck Collaboration XV 2013) to derive parameter con- 
straints in Planck Collaboration XVI (2013). In Sect. 6.1 we 
highlight a subset of parameter results where the information 
provided by the lensing likelihood has proven particularly use- 
ful. In the concordance ACDM cosmology, there is believed to 
be a correlation between the CMB lensing potential and the low- 
i temperature anisotropics, driven by the eff'ects of dark energy. 
We also present a measurement of this correlation in Sect. 6.2. 
Finally, we connect our lensing potential map to other tracers of 
large-scale structure with several illustrative cross-correlations 
using galaxy, quasar, cluster and infrared source catalogues in 
Sect. 6.3. These main results are followed in Sect. 7 by a large 
suite of systematic and consistency tests, where we perform null 
tests against a variety of diff'erent data cuts and processing. We 
conclude in Sect. 8. A series of appendices provide further de- 
tails on some technical aspects of our methodology and lensing 
potential estimates. 

Throughout this paper, when we refer to the concordance or 
fiducial ACDM cosmology we are referring to a model with 
baryon density ^^^b = ^hh^ = 0.0480, cold dark matter den- 
sity ojc = ^ch^ = 0.1199, density parameter for the cos- 
mological constant ^Ia = 0.6910, Hubble parameter Hq = 
lOOh km s"^ Mpc"^ with h = 0.6778, spectral index of the power 
spectrum of the primordial curvature perturbation n^ = 0.96, am- 
plitude of the primordial power spectrum (at k = 0.05 Mpc"^) 
As = 2.21 X 10"^, and Thomson optical depth through reion- 
ization r = 0.093. These values were determined from a pre- 
publication analysis of the Planck temperature power spectrum, 
but are consistent with the best-fit values quoted in Planck 
Collaboration XVI (2013). 



2. Methodology 

In this section, we detail our methodology for reconstructing 
the lensing potential and estimating its angular power spectrum. 
These are both accomplished by exploiting the distinctive statis- 
tical properties of the lensed CMB. 

(I) If we consider a fixed lensing potential applied to multi- 
ple realizations of the CMB temperature anisotropics, then 
lensing introduces statistical anisotropy into the observed 
CMB; the fluctuations are still Gaussian, however the co- 
variance varies as a function of position and orientation 



on the sky. We use this idea to obtain a (noisy) estimate of 
(p(h). The noise of this map is a combination of instrumen- 
tal noise and statistical noise due to the fact that we only 
have a single realization of the CMB to observe, analo- 
gous to shape noise in galaxy lensing. 
(II) If we consider averaging over realizations of both the lens- 
ing potential and the CMB fluctuations, then lensing intro- 
duces non-Gaussianity into the observed CMB. This ap- 
pears at lowest order in the connected part of the CMB 
4-point function, or trispectrum^ . We use this to measure 
the lensing power spectrum C^*^. 

The estimators that we use are derived from maximizing the like- 
lihood function of the lensed CMB under the hypotheses above, 
and should be optimal (in the minimum- variance sense). In cases 
where we have made suboptimal choices, we provide estimates 
of the loss of signal-to-noise. 

2.1. Lens reconstruction 

To gain intuition for the process of lens reconstruction, it is use- 
ful to consider the efl'ect of lensing on a small patch of the sky. 
Lensing remaps the temperature fluctuations by a deflection field 
V(p(h). The part of V(p(h) that is constant over our patch is not 
an observable eff'ect; it describes only a re-centering of the map. 
The variation of the deflection field across the patch is observ- 
able, however. This can be usefully decomposed into conver- 
gence (k) and shear modes (7+, y.) as 



■ V,-V,-0(/z) : 



K-\-y+ 7_ 
7- K-y+ 



in). 



(4) 



If we observe a patch that is small enough that these quan- 
tities can be taken as constant, then the observational con- 
sequences are simple. The convergence mode causes a local 
change of scale, either magnifying or demagnifying the fluctu- 
ations. Taking the local power spectrum of our small patch, we 
would find that the CMB peaks would shift to larger or smaller 
scales, relative to the full- sky average. The shear modes also de- 
scribe changes of scale, however they are now orientation depen- 
dent. On a small patch, convergence and shear estimators can be 
constructed from local estimates of the (orientation-dependent) 
power spectrum and then stitched together to recover the lensing 
potential (Zaldarriaga & Seljak 1998; Bucher et al. 2012). This 
procedure describes a quadratic estimator for the local conver- 
gence and shear. 

From the description above, it is not immediately clear how 
to go about stitching together estimates of convergence and shear 
in diff'erent regions of the sky, or what weight to give the local 
power spectrum estimates as a function of scale. These questions 
can be resolved by considering a generic form for the quadratic 
estimator, and optimizing its weight function for sensitivity to 
lensing (Okamoto & Hu 2003). To first order in the lensing po- 
tential, the statistical anisotropy introduced by lensing appears 
as an off'-diagonal contribution to the covariance matrix of the 
CMB: 



LM ^1 mi ,^2^2 ^ ' 



(5) 



^ The ISW-lensing correlation also introduces a non-zero bispectrum. 
When correlating the reconstructed 0(n) with the large-angle tempera- 
ture anisotropics in Sect. 6.2, we are probing this bispectrum. 
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where the average () is taken over CMB realizations with a fixed 
lensing potential. Here the bracketed term is a Wigner 3 j sym- 
bol, = / d^nYlj^(h)(p(h) is the harmonic transform of the 
lensing potential, and the weight function W^^^^l given by 



Here Cj^ is the theoretical power spectrum of the lensed CMB. 
Note that we use the lensed power spectrum here, rather than 
the unlensed spectrum that is sometimes used in the literature, 
as this is accurate to higher order in (Lewis et al. 2011), an 
improvement which is necessary at Planck sensitivity (Hanson 
et al. 2011). 

Now we construct a quadratic estimator to search for the co- 
variance which is introduced by lensing. We will use several dif- 
ferent estimators for the lensing potential, as well as to probe 
possible point source contamination, and so it will be useful to 
keep this discussion as general as possible. A completely generic 
quadratic estimator for the lensing potential can be written as 

^- = EMi„,,[w-x-], (7) 

L'M' 

where is a normalization matrix, and xlm is a quadratic 
"building block" which takes in a pair of filtered sky maps ff^ 
and ff^, and sums over their empirical co variance matrix with a 
weight function ^^^^^l- 



2 Zj ^ ^ \mi m2 -M ] 



T^}\ T?l . (8) 



The "mean-field" term x^^ accounts for all known sources of 
statistical anisotropy in the map, which could otherwise bias the 
lensing estimate. It is given by 



-MF 



- i V (-1)^1^^ ^2 L \ -(1) 



tZ^^ (9) 



where the ensemble average here is taken over realizations of the 
CMB and noise. 

We may now optimize the generic quadratic estimator above. 
If the primordial CMB fluctuations and instrumental noise are 
Gaussian and the lensing potential is fixed, then the likelihood 
for the observed CMB fluctuations is still a Gaussian, which may 
be maximized with respect to the lensing potential modes (piM 
(Hirata & Seljak 2003a). The optimal quadratic estimator is the 
first step of an iterative maximization of this likelihood, and it 
has been shown that additional iterations of the estimator are not 
necessary for temperature lens reconstruction (Hirata & Seljak 
2003a; Okamoto & Hu 2003). The optimal quadratic estimator 
has the following choices for the weight function and filtering. 

(I) The weight function should be a matched filter for the 
covariance induced by lensing (i.e., one should use 0, with 
weight function given by Eq. 6). We shall use this weight 
function for all of our fiducial results, although for consis- 
tency tests we will also use "bias-hardened" estimators, 
which have weight functions constructed to be orthogo- 
nal to certain systematic eff'ects (Namikawa et al. 2012a). 
This is discussed further in Sect. 7.4. 



(II) The filtered temperature multipoles T^m should be given 
by T{fn = (C~^T){ni, where T is a beam-deconvolved sky 
map and C is its total signal-hnoise covariance matrix. We 
describe our approximate implementation of this filtering 
in Appendix B. When combining multiple frequencies for 
our minimum- variance estimator, all of the available data 
is combined into a single map which is then filtered and 
used for both input multipoles of the quadratic estimator. 
It can be desirable to use diff'erent pairs of maps however, 
and we use this for several consistency tests. For exam- 
ple, we feed maps with independent noise realizations into 
the quadratic estimator to avoid possible noise biases in 
Sect. 7.3. 

In the quadratic maximum-likelihood estimator, the mean-field 
correction emerges from the determinant term in the likelihood 
function, and it can be seen that the normalization matrix *R is 
the Fisher matrix for the 0lm; this means that the normalization 
is the same as the covariance matrix of the lens reconstruction, 
and so the unnormalized lensing estimate = x - x^^ is equiv- 
alent to an inverse- variance- weighted lens reconstruction, which 
is precisely the quantity needed for most statistical analysis. This 
is why we have denoted it with an overbar, in analogy to T. 

We choose to treat the map noise as if it were homogeneous 
when constructing the filtered f^^, and do not account for vari- 
ation with hit count across the sky. This is a slightly suboptimal 
filtering choice; in Appendix. B we estimate that it leads to a 
5% loss of total signal-to-noise when constraining the power 
spectrum of the lensing potential. The advantage of this ap- 
proach, however, is that far from the mask boundaries our fil- 
tering asymptotes to a simple form, given by 

Tim ~ [C'J^ + C^^j T^jn = F{T{fn, (10) 

where Cj^ is the temperature power spectrum and C^^ is the 
power spectrum of the homogeneous noise level that we use in 
our filtering. For the purposes of compact notation, in the fol- 
lowing equations we combine both of these elements in the "fil- 
ter function" F^. The asymptotic form of our filtering, Eq. (10), 
will prove useful, as it means that the normalization of our esti- 
mator, as well as its variance and response to various systematic 
eff'ects, may be accurately modelled analytically. It allows us to 
propagate uncertainties in the beam transfer function and CMB 
power spectrum, for example, directly to our lens reconstruction. 
This filtering choice also means that the normalization does not 
vary as a function of position on the sky, which simplifies the 
analysis of cross-correlations between the lensing potential map 
and external tracers. Under the approximation of Eq. (10), the 
estimator normalization is given by 



(11) 



where the response function for filtered maps T^^^ and T^^^ 



IS 



^x0,(l)(2) 



(2L+l)f-^2 



T.\KeALp'i:p'^- (12) 



This can be read as "the response of estimator x to lensing on 
scale L". The filter functions are those used for f and f 
respectively. In cases where the filter functions are obvious, we 
will drop the indices above. 

Putting all of the above together, for a chosen quadratic esti- 
mator X we obtain normalized, mean-field-debiased estimates of 
the lensing potential as 



(13) 
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Note that our normalization function is only approximate, 
but we will verify its accuracy in Sect. 4. For the standard lens- 
ing estimator of Okamoto & Hu (2003) (which uses the weight 
function of Eq. 6), we use x = (p. This estimator is denoted sim- 
ply as 0LM. 



2.2. Lensing power spectrum estimation 

We form estimates for the power spectrum of the lensing po- 
tential by taking spectra of the lensing estimates from Sect. 2.1, 
using a simple pseudo-Q estimator. In order to reduce mode 
coupling, as well as to downweight regions near the analysis 
boundary where the mean-field due to masking can be large, we 
take the power spectrum from an apodized version of our lensing 
estimate, given by 



L'M' 



YL^M'(n)PL'(p] 



L'M' 



, (14) 



where M(n) is an apodized version of the analysis mask M(h) 
used in our filtering and Pl = L(L -h 1) is an approximate pre- 
whitening operation. The construction of M(fi) is described in 
Sect. 3. Our fiducial apodization occurs over a band of approxi- 
mately 5°, and eff'ectively reduces the sky fraction by 9%. 

The power spectrum of probes the 4-point function of 
the observed CMB, which contains both disconnected and con- 
nected parts. We model it as being due to a combination 
of Gaussian CMB fluctuations, lensing efl'ects and unresolved 
point-source shot noise, and estimate the power spectrum of the 
lensing potential with 



^L,x ~ 



L INO 



M 



AC 



,001 

L InI 



AC 



4(f) \ 

L IpS 



AC 



4(f) \ 

L ImC ' 



(15) 



where /sky,2 = / dhM^(h)/4n is the average value of the square 
of the apodizing mask. The first line of Eq. (15) isolates the con- 
nected part of the CMB 4-point function, or trispectrum, which 
would be zero for Gaussian fluctuations. The second line con- 
tains corrections which isolate the part of the trispectrum which 
is directly proportional to the non-Gaussianity induced by C^^. 
In the following paragraphs, we explain these terms in more de- 
tail. 

The first correction term AC^'^Ij^q subtracts the (large) dis- 
connected contribution to the power spectrum of 0. To determine 
this term, we use the data-dependent subtraction which emerges 
for maximum-likelihood estimators of the CMB trispectrum 
(Regan et al. 2010; see also Appendix D). For lensing, this pro- 
cedure has the additional advantage of reducing the correlation 
between diff'erent multipoles L ^ U of the lens reconstruction 
(Hanson et al. 2011), as well as reducing sensitivity to uncer- 
tainties in our model of the CMB and noise covariance matrices 
(Namikawa et al. 2012a). Writing the power spectrum of cpiM 
explicitly as a function of the four inverse- variance filtered tem- 
perature maps 



^# |- j-(l) j-(2) j-(3) j-(4)j ^ 



2L 



(16) 



the disconnected contribution reads 



AC 



l,xIno 



^L,X MC ^ MC" ^ MC" ^ MCJ 



^# r^(l) ^(2) f(3) f(4)],^(f>(f> \f(l) f(2) f(3) f(4)] 

^L,jc[ MC'^ '^MC'^ J ^L,X [ MC ^ ' ^ McJ 

_,r^\f(l) f'(2) f(3) f'(4)] , ff-Cl) f(2) f(3) f(4)] 

^L,x[ '^MC'^ ' ^ MCj ^L,X [ '^MC'^MC'^ J 



_ ^# r^(l) f(2) f(3) f(4) ] 

^L,X MC ^ MC" ^ MC ^ MC'J 



, (17) 



where Tmc indicates a Monte-Carlo simulation of the corre- 
sponding map. The ensemble average is taken over two sets of 
independent realizations MC and MC'. Note that because of the 
way we have used pairs of Monte-Carlo simulations and data 
with independent CMB and noise realizations, the mean-field 
correction is zero for all of the terms above. 

The term AC^'^Ini corrects for the 'W^^^" bias due to sec- 
ondary contractions of the lensing trispectrum (Hu 2001; Kesden 
et al. 2003). It is only a large eff'ect at L > 100, and so we calcu- 
late it using a flat-sky expression as 



Ac; 



1 



1^1,^1 



Ni aDX(p,(m)qDX(P,(3)(4) 



r (fh r £i 

J (2^ 



(l)^(2)^(3)^(4)^x 

<^^fi>^(-^l' ^3)^^-/2, /4) 

+ ^^f^'|'^^(-^l'^4)W^(-/2,/3) 



(18) 



where h + h = h + h = L and C 



'(f)(f), fid. 



is a fiducial model for 



the lensing potential power spectrum. The W(l, V) are flat- sky 
analogues of the full-sky weight functions. The flat- sky lensing 
weight function, for example, is 



L + C 



TTj 

\l2\ ^2 



(19) 



M 



The N^^"* term is proportional to the lensing potential power 
spectrum, and so in principle it should be used to improve our 
constraints on C^^ rather than subtracted as an additive bias. 
However the statistical power of this term is relatively small 
at Planck noise levels. From a Fisher matrix calculation, the 
trispectrum contractions which source the A/^^^^ term are only de- 
tectable in the Planck data at 4cr significance, compared to the 
approximately 25cr for the primary contractions. We choose sim- 
ply to subtract the A/^^^^ term from our power spectrum estimates. 
There is a small cosmological uncertainty in the A/^^^^ correction 
due to uncertainty in the C^^ power spectrum, which we discuss 
in Sect. 5.3. 

The AC^'^lpg term is a correction for the bias induced by the 
non-Gaussianity of unresolved point sources, discussed in more 
detail in Sect. 2.4. 

Finally, the AC^'^lj^^ term is a small correction that we ob- 
tain by estimating C^^ following the procedure above on a num- 
ber of lensed CMB realizations, and then subtracting the input 
power spectrum. This term can be non-zero due to pseudo-Q 
leakage eff'ects from masking, which we have not accounted for 
other than apodization, errors in our calculation of the A/^^^^ term, 
or errors in the normalization at the power spectrum level. We 
will find that AC^'^lj^^ is sufficiently small that in practice it 
does not matter whether we account for it as a renormalization 
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or an additive offset, and we choose to treat it as an offset for 
simplicity. 

We will ultimately characterize the uncertainty of C^j^ by 
Monte Carlo, however the following analytical expression is a 
useful approximation 

y^Kci'j " = ^ ^ [^L + <]' (20) 



where A^^'^^ is the reconstruction noise level. We take it to be 



^'"^ (2L + 1) <^^0,(l)(2)<^x0,(3)(4) 



2^rr,(i)(3)^rr,(2)(4)^ (21) 



where C^'^^^^ is the ensemble-average cross-spectrum between 



r^^^ and T^^\ given by 



C 



TT,(i)U) 



2L 



2.4. Unresolved point-source correction 

At the high multipoles (£ > 1000) that provide most of the modes 
which are useful for lens reconstruction, the contribution from 
unresolved (and therefore undetected and unmasked) extragalac- 
tic foregrounds becomes apparent in the Planck power spectra. 
The simulation and analytical modeling work of Osborne et al. 
(2013) indicates that at Planck resolution and sensitivity, the 
contribution of unresolved foregrounds to lens reconstruction is 
small but potentially non-negligible, with the largest contribu- 
tion coming from the shot noise trispectrum of unresolved point 
sources. 

Our approach to this bias is to measure the amplitude of the 
shot noise contribution and correct the measured lensing spec- 
trum. We measure the shot noise amplitude using the power 
spectrum of a quadratic estimator designed to detect the "noise" 
due to sources, similar to the approach advocated in Munshi et al. 
(2011a). The trispectrum (or connected four point function) for 
point source shot noise is defined in position space as 



77Zi(^iM^^*>' (22) {T{h,)T{h2}T{h,)T{h,)) 

M rf4c/. 



where fskyXOU) fraction of sky common to both T^^^ and 

T^^\ We will use this analytical estimate of the variance for 



weighting our C^^ estimates in statistical analyses. 



2.3. Lensing power spectrum likelihood 

Based on our measurements of the lensing potential power spec- 
trum, we construct a Gaussian likelihood based on bins in C^^. 
Our likelihood has the form 

- 2 In £^(C^^) = B\ (Cf - ) [s-i]'' (C^^ - C^^) , 

(23) 

where ^ represents a binning function, E is the covariance matrix 
between bins, and sums are performed over paired upper/lower 
indices. As the shape and amplitude of C^^ are constrained 
strongly in the concordance ACDM model, with which our re- 
sults are broadly consistent, we choose a binning function de- 
signed to maximize our sensitivity to small departures from the 
fiducial ACDM expectation. This is given by 



,fid. 



ifL^^<L<L^,„ 



(24) 



where Vl is defined in Eq. (20) and L|^.^ and V^^^ define the 
multipole range of the bin. Our binned results correspond to 
an estimate of the amplitude of a fiducial C^^ power spectrum 
in a given multipole range, normalized to unity for the fidu- 
cial model. We denote these amplitude estimates explicitly for 
a given lens reconstruction Cf^ as 



(25) 



In principle there are several reasons why a likelihood ap- 
proach such as the one we have described above could fail; 
the usual issues with a pseudo-Q likelihood, such as the non- 
Gaussianity of the C^^ are compounded by the fact that the 
estimates themselves are non-Gaussian, and derived from the 
temperature data itself, which means that the measurement un- 
certainties on the lens reconstruction are potentially correlated 
with those of the temperature power spectrum. In Appendix D, 
we validate the above approach to the likelihood by considering 
these issues in more detail. 



S 6(h\ - h2)S(h2 - h3)S(hi, - W4), (26) 
or in harmonic space: 

(T^^mi T£^m2 T^hm3 T^^m^ ^ ~ dflY^^mi ^^^2^2 ^hm'i ^^^4 • (27) 

The quadratic estimator that we use to measure the amplitude S ^ 
is denoted as s, and is defined by the weight function 



(2A + l)(2^2 + l)(2L+l)/^i £2 L 



'^i^2L - ^^ 4^ \^0 0/' 

We measure the amplitude of the shot-noise trispectrum 5" ^ as 



(28) 



-1 



(29) 



where C^^ is defined analogously to Eq. (15) and the sums are 
taken over L,L 
S lis given by 



taken over L,U e [L^.^,L^^^]. The power spectrum weighting 



Sl = 



1 



1 



/sky,2 2L + 1 



.^^,(1)(2)<^^^,(3)(4) • 



(30) 



In the limit that [^^^j^, ^max] - [0' ^1 this estimator is equiva- 
lent to measuring the shot noise trispectrum using the 1 -point 
kurtosis of the map. The advantage of the trispectrum-related 
power spectrum approach is that it allows us to separate out re- 
gions of the trispectrum measurement which are contaminated 
by lensing, as well as to look for evidence of clustering, which 
would appear as a deviation of the measured spectrum C^^ from 
the shape expected for unclustered sources. We then calculate 
and remove an estimated bias to the measured lensing potential 
given by 



3. Data and cuts 

Planck sky maps The majority of the results in this paper 
are based on the Planck nominal mission frequency maps at 
100, 143, and 217 GHz, built from the first 15.5 months of 
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data. These are in HEALPix^ format, with resolution param- 
eter A^side = 2048, corresponding to pixels with a typical 
width of V.l. These have effective noise levels of approxi- 
mately lOSyuKarcmin at 100 GHz, 45yuKarcmin at 143 GHz, 
and 60yuKarcmin at 216 GHz. The beam widths shrink with 
frequency, and are 10' at 100 GHz, 1' at 143 GHz, and 5' at 
217 GHz. Our primary products - a map of the lensing poten- 
tial and an associated power spectrum likelihood - are based 
on a minimum- variance (hereafter MV) combination of the 143 
and 217 GHz maps. Although lensing may be detected at a sig- 
nificance of approximately lOcr at 100 GHz, the CMB modes 
used in this reconstruction are already sample-variance domi- 
nated at 143 and 217 GHz, and so adding it to the MV com- 
bination would lead to negligible improvement. In addition to 
these maps, we also use the 857 GHz Planck map as a dust 
template, which is projected out of both maps independently 
in our filtering procedure (described further in Appendix B). 
This projection is primarily intended to remove diff'use Galactic 
dust contamination, although it also removes a portion of the 
GIB fluctuations, which have a similar spectral index to that of 
Galactic dust over these frequency bands Planck Collaboration 
XXIV (2011); Planck CoUaboration XVIII (2011). As a sim- 
ple approach, this template projection overlooks several poten- 
tial difl&culties, including variation of Galactic dust spectral in- 
dices across the sky as well as mismatch between the beams at 
100/143/217 and 857 GHz, although we find it is adequate for 
our purposes. As lens reconstruction is most sensitive to small- 
scale modes, the coupling to large-scale Galactic foregrounds 
is relatively weak. In Sect. 7 we will also perform lens recon- 
struction using the more rigorously component-separated maps 
(Planck Collaboration XII 2013), finding good agreement with 
our baseline results. 

In analyzing these maps, we use the fiducial beam trans- 
fer functions described in Planck Collaboration VII (2013) and 
Planck Collaboration ES (2013). There are uncertainties asso- 
ciated with these transfer functions, which we propagate to an 
uncertainty in the normalization of lensing estimates in Sect. 5.2. 

Galaxy mask We avoid the majority of Galactic foreground 
power using the temperature analysis masks described in Planck 
Collaboration ES (2013). These are constructed using a combi- 
nation of the 30 GHz and 353 GHz maps, corrected for a CMB 
estimate and smoothed to 5' and thresholded until a desired sky 
fraction is obtained (Planck Collaboration XII 2013). For our 
baseline results, we use the 70% masks (which remove 30% of 
the sky), although in Sect. 7 we will show that we obtain con- 
sistent results with both larger and smaller masks. When com- 
puting the power spectrum, we additionally multiply our es- 
timates by an apodized version of the Galaxy mask. Each pixel 
out of the masked region, is multiplied by a weight varying be- 
tween zero at the mask boundary and unity at a distance greater 
than 5° from the closest masked pixel. We use a sinusoid weight 
function, similar to the one used in Namikawa et al. (2012a) and 
Benoit-Levy et al. (2013). Note that due to this apodization, the 
eff'ective sky fraction used in our power spectrum analysis is ap- 
proximately 9% lower than the sky fraction of our reconstruc- 
tion. 

CO and extended-object masks We mask regions believed to 
be contaminated by carbon-monoxide (CO) emission lines at 
100 and 217 GHz using the type 2 CO map described in Planck 

^ http://healpix.jpl.nasa.gov 



Collaboration XIII (2013). We threshold all pixels above 60 yuK 
in the map after degrading to /"/side = 256, and discard by hand all 
isolated pixels which are removed by this cut but do not appear to 
be extended CO regions. This low-resolution mask is prograded 
back to A^side = 2048 and smoothed at 20', slightly larger than the 
map resolution of 15'. We additionally exclude extended nearby 
objects (the two Magellanic Clouds, M31, M33, and M81) by 
cutting out disks centred on their locations. The radii of these 
disks range from 250' for the Large Magellanic Cloud to 30' for 
M81. 

Point-source masks We remove detected point-sources (or oth- 
erwise compact objects) using a mask constructed from a combi- 
nation of sources identified in the Planck Early Release Compact 
Source Catalogue (ERCSC; Planck Collaboration Int. VII 2013), 
the SZ clusters from the Planck Cluster Catalogue (PCC; 
Planck Collaboration XXIX 2013), and the Planck Catalogue of 
Compact Sources (PCCS; Planck Collaboration XXVIII 2013), 
We produce individual masks for 100, 143, and 217 GHz, cutting 
out detected sources with disks having radii of either 3cr or 5cr 
depending on their flux level, where a describes the Gaussian 
beam- width of the given channel and is taken to be 4.1' at 
100 GHz, 3.1' at 143 GHz, and 2.1' at 217 GHz. For the PCCS 
we mask all sources with SNR > 5 in a given band, or SNR > 10 
in either adjacent frequency band, using a 3(T disc, extending to 
a 5(T disc for any sources with SNR > 10 in the target band. For 
the ERCSC, we make a 3(T cut for every source detected at the 
target frequency, as well as for every source which is detected 
with a flux greater than 1 Jy at 70 GHz, 0.7 Jy at 100 GHz, 0.4 Jy 
at 143 GHz, 0.5 Jy at 217 GHz, or 0.9 Jy at 353 GHz. For sources 
that satisfy the brightness criterion at the target frequency, we 
also enlarge the cut to 5cr. For SZ clusters we make a cut be- 
tween 3(T and 5cr depending on the cluster angular size ^500. We 
further mask the cool cores (CC) using the ECC catalog in a 
similar fashion, using a maximum 15' excluding radius. 

In Sect. 7.2 we will test our sensitivity to this SNR cut, vary- 
ing it from SNR > 5 to SNR > 4.2 and SNR > 10. 



4. Simulations 

We require simulations of the data described in Sect. 3 for sev- 
eral aspects of our analysis: to determine the lensing-mean field 
0^^, to determine the AC^^^I^^^ and AC^'^lj^^ correction terms 
of Eq. (15), and to determine our measurement error bars, as well 
as to validate our reconstruction methodology. Here we outline 
the basic algorithm we use to generate these simulations. 

(1) Simulate the unlensed CMB and lensing potential. Using the 
fiducial ACDM theoretical unlensed power spectra for the 
cosmology given on Pg. 4, we generate realizations T^^ and 
(pLM of the unlensed temperature and lensing potential with 
appropriate correlations. By default we simulate these fields 
up to ^max = 2560. In addition to these stochastic compo- 
nents, we also include a dipolar lensing mode in our sim- 
ulations, to account for the Doppler aberration due to our 
motion with respect to the CMB (Kamionkowski & Knox 
2003; Challinor & van Leeuwen 2002): 

(Pp{h)^P-h, (32) 

where P is in the direction of Galactic coordinates (/, b) = 
(264.4,48.4), with amplitude \/3\ = 0.00123, corresponding 
to a velocity of 369 km s"^(Hinshaw et al. 2009). 
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(2) Lens the temperature field. We use a fast spherical harmonic 
transform to compute the temperature and deflection fields 



T'"'\ii) = Y,YumT2^ 

im 

d{h) = Yj^YLM{h). 



(33) 



This computation is performed on an equicylindrical pix- 
elization (ECP) grid with {Nq^N^) = (16384,32768) 
equally- spaced points in (6,(p). We determine the lensed 
temperature on this ECP grid as 



T(h) = T''''\h + d{n)). 



(34) 



This equation is evaluated following the same algorithm as 
the LENSPix code (Lewis 2005), with 2D cubic Lagrange 
interpolation on the ECP grid, which has conveniently uni- 
form pixel spacing. We then obtain the lensed temperature 
multipoles T^^ by a reverse harmonic transform. The aver- 
age power spectrum of these lensed simulations agrees with 
the implementation in CAMB (Challinor & Lewis 2005) to 
better than 0.1% at / < 2048. We have also used the algo- 
rithm of Basak et al. (2009), which gives consistent results. 
In addition to the lensing dipole, our motion with respect 
to the CMB also causes a small modulation of the observed 
fluctuations, which we determine from the lensed tempera- 
ture field as 

T^'\h) = fyfi ' hT(h), (35) 

where fy is a frequency dependent boost factor that we ap- 
proximate as 1.5, 2, and 3 at 100, 143, and 217 GHz, re- 
spectively. There is an extended discussion of this factor in 
Planck CoUaboration XXVII (2013). 
(3) Generate a Gaussian realization of extragalactic foreground 
power r^^. We simulate Gaussian foreground power at 100, 
143, and 217 GHz such that the auto- and cross-spectra of 
the foreground power at these frequencies has the form 



r^xB _ 4PS ^PS , jCIB ^CIB 



(36) 



where A and B label frequency bands. The template power 
spectra are based on the modelling of Planck Collaboration 
XV (2013); Planck CoUaboration XVI (2013), and given by 



P3 27t 
' "30002' 



^CIB 



In 
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We use the parameters 
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PS 
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PS 
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= 43yuK2 



PS 

100x217 
PS 

217x217 
CIB 



(37) 



= 44yUK^ 
= 57yUK^ 



= 14a/K^ A--2i7 = 54yuK 



The parameters above were obtained by fitting the auto- 
and cross-spectra of our inverse-variance-filtered tempera- 
ture maps to the templates above. Our point source masks 
difl'er from those used in Planck Collaboration XVI (2013), 
and so these numbers cannot be directly compared to the val- 
ues there. The purpose of including a Gaussian foreground 
component in our simulations is primarily so that our mea- 
surement error bars will include the additional scatter due to 
the foreground power (we will see in Sect. 5 that foreground 
power constitutes a few percent of our total error budget). 
Our lensing power spectrum results are otherwise unaff'ected 
by this component of our simulations. 



(4) Convolve j'/'+^fg +7^0(1 ^j^^ (asymmetric) instrumental 
beam using Eq. (C.7). This results in a HEALPix map of the 
beam-convolved CMB T^(h), from which we extract beam- 
convolved multipoles T^^. 

(5) Make a final map of the observed sky T(h) with the 
HEALPix A^side = 2048 pixelization, with 



T(h) = Y^m 



(38) 



where B^^^ is the HEALPix A^side = 2048 pixel window 
function, which accounts for the eff'ective smearing of the 
beam due to the binning of hits into pixels. Although this 
smearing is an inadequate description of the pixelization- 
smearing eff'ect for any individual pixel, it is a good approx- 
imation when averaged over all pixels on the sky. We will 
discuss pixelization eff'ects further in Sect. C.4. The quantity 
n(h) is a simulated realization of the instrumental noise. We 
use noise realizations from the FFP6 simulation set (Planck 
Collaboration I 2013). In order to make the power spectrum 
of our simulations more closely match the power spectrum 
of the data, we find that it is necessary to include an ad- 
ditional white noise component with an amplitude of 20, 
10 and lOyuKarcmin at 100, 143 and 217 GHz respectively. 
Indeed the FFP6 simulation are produced on the basis of the 
half-ring estimates of the noise that are known to underesti- 
mate it by about 2% (Planck Collaboration VI 2013) . 

The simulations described above are used throughout this work 
both to debias our lensing estimates, as well as to characterize 
their uncertainty. Here we briefly use them to validate the nor- 
malization of our lensing estimates. There are two relevant nor- 
malizations: the map-level normalization which describes how 
our reconstruction traces the input realization, as well as the 
analogous normalization for the lensing potential power spec- 
trum. Throughout this work, we use the analytical normalization 
estimates of Eq. (12) to normalize our map estimates. In Fig. 3 
we plot the average cross-correlation between the input and re- 
constructed (pLM for our minimum-variance reconstruction, di- 
vided by a simple /sky factor to account for missing power in the 
mask. We can see that our analytical map-level normalization is 
accurate at the 1% level. This is small enough that we do not 
further renormalize our map. At the power spectrum level, we 
find similar accuracy of our normalization in the 40 < L < 400 
region which forms the basis of our fiducial lensing likelihood. 
The power spectrum estimates in Fig. 3 are obtained by averag- 
ing the first line of Eq. (15) over a set of lensed but otherwise 
Gaussian simulations. This 1% error in the power spectrum nor- 
malization is accounted for in our likelihood by the AC^'^lj^^ 

term. Given the good agreement of the C^^ power spectrum 
used in our simulations with that measured from data, it does 
not matter whether we incorporate the Monte-Carlo correction 
as a multiplicative or additive term. We have chosen to incorpo- 
rate it as additive for simplicity. At higher multipoles, there is 
additional scatter in our power spectrum normalization test (up 
to 10% for the plotted bins). It is difficult to know from this test 
whether it is due to actual scatter in our power spectrum normal- 
ization (which seems unlikely, given that the map-level normal- 
ization is quite accurate), a small misestimation of the N^^'^ term 
(which begins to dominate over C^^ above L = 500), or simply 
pseudo-Q leakage issues. In any case, it is still small, and is ab- 
sorbed into our AClImc t^rm when estimating power spectra at 
these multipoles. Although we have only presented results for 
our minimum- variance reconstruction here, those for individual 
frequency results are similar. 
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As a visual illustration, and a preview of our data results in 
Sect. 6, in Fig. 4, we show a simulated lens reconstruction as 
well as the input map, which gives a visual impression of the 
signal-to-noise in our lens reconstruction. 
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Fig. 3. Validation of our estimator normalization for simulations 
of the MV reconstruction at the map and power spectrum levels. 
The map normalization (plotted as 0^^^) is tested by taking the 
cross- spectrum of the input with the reconstruction averaged 
over Monte Carlo simulations, divided by an /sky factor to ac- 
count for missing power in the mask. The power spectrum nor- 
malization (plotted as 00) is obtained by averaging the first line 
of Eq. (15) over simulations, and then comparing it to the ex- 
pected value, which is C^j^ + /S.C^j^\^^ because our simulations 
do not contain point source non-Gaussianity. 
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Fig. 4. Simulation of the Wiener- filtered lensing potential esti- 



mate 



lWF - 



^ C7(cPlm 



IMF 



) for the MV reconstruction (left). 



and the input realization (right; filtered by C^^H^^ to be di- 
rectly comparable to the Wiener estimate). Both maps show the 
southern Galactic sky in orthographic projection. The lensing 
reconstruction on the data is noise dominated on all scales, how- 
ever correlations between the two maps can still be seen visually. 



5. Error budget 

In this section, we describe the measurement and systematic er- 
ror budget for our estimation of the lensing potential power spec- 
trum. This is broken down into three sections; in Sect. 5.1 we 
describe our measurement (or "statistical") error bars, which are 
due to the fact that we have only a single noisy sky with a finite 
number of modes to observe. In Sect. 5.2 we consider uncer- 
tainty in the instrumental beam transfer function, which we will 
see propagates to a normalization uncertainty for our lensing es- 
timates. Finally, in Sect. 5.3 we discuss the eff'ect of cosmolog- 
ical uncertainty; possible errors in the fiducial model for 
result in a normalization uncertainty for our lensing estimates, 
and uncertainties in the fiducial C^j^ power spectrum lead to un- 
certainties in the A^^^^ correction. As a guide to the relative size 
and scale dependence of these terms, in Fig. 5 we summarize the 
error budget for our fiducial minimum- variance lens reconstruc- 
tion, based on 143 and 217 GHz. Individual frequency bands, as 
well as 100 GHz are qualitatively similar. 



5.1. Measurement 

Although our measurement uncertainties are ultimately assigned 
by Monte Carlo, we can use the analytical expression of Eq. (21) 
to gain intuition for how they are sourced by various compo- 
nents. Our simple model of the sky after masking and dust clean- 
ing is that it consists of three uncorrelated signals: CMB, instru- 
mental noise, and unresolved foreground power. The noise vari- 
ance of the lens reconstruction in Eq. (21) involves two power 
spectra, and so we can think of the noise contribution as the 
sum of six possible terms involving pairs of the CMB, noise, and 
foreground power spectra. In Fig. 6 we combine these contribu- 
tions into three representative contributions to the reconstruction 
noise: "pure CMB" in which both spectra are due to CMB fluc- 
tuations; the "noise" contribution in which either both spectra 
are those for noise power, or one is noise and one is CMB; and, 
finally, the "foreground" contribution in which either one or both 
of the spectra are due to unresolved foreground power. We can 
see that for most reconstruction multipoles, the pure CMB con- 
tribution constitutes the largest part of the reconstruction noise, 
followed by noise. The unresolved foreground power is a fairly 
small contribution to our measurement error. Note that the dom- 
inant terms for both the "noise" and "foreground" contributions 
are the ones in which one of the spectra is a CMB fluctuation. 
For this reason, we will focus less on the use of cross- spectra to 
avoid noise biases than is done for the usual CMB power spec- 
tra (Planck Collaboration XV 2013), although we will perform 
consistency tests using cross- spectra of data to avoid noise bi- 
ases. Note that our realization-dependent method for removing 
the disconnected noise bias (Eq. 17) means that the majority of 
this contribution is estimated directly from the data itself, re- 
ducing our sensitivity to uncertainty in the noise and foreground 
power. 



5.2. Beam transfer function 

Errors in the efl'ective beam transfer function appear as an error 
in the normalization of our lensing estimates. For simplicity here 
we will describe the case for a single standard quadratic lensing 
estimator that uses the same map for both of its inputs, although 
when dealing with combinations of channels for our actual re- 
sults we account for difl'erences in the beam transfer function 
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Fig. 5. Relative contribution of statistical measurement uncer- 
tainty (grey band, representing the Icr uncertainty on C^j^ for the 
plotted bins) and the systematic errors we assign due to uncer- 
tainty in the N^^"^ bias correction, point-source shot noise correc- 
tion, and the beam transfer function. The contributions represent 
error eigenmodes in C^^ , and so are completely correlated be- 
tween L. Note that these contributions are much smaller than 
our measurement uncertainty, and they have been multiplied by 
a factor of 10 for clarity. 
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Fig. 6. Relative contributions of CMB, instrumental noise, and 
foreground power terms discussed in Sect 5.1 to the approxi- 
mate lens reconstruction variance of Eq. (21), previously plotted 
in Fig. 1. As discussed in Sect. 5.1, the finite number of CMB 
modes which we observe is the dominant source of variance for 
the lens reconstruction. 



and errors between "legs" of the estimator (or pair of estimators, 
in the case of the lensing potential power spectrum). 

We model the eff'ect on from a discrepancy between the 
"measured" and "true" eff'ective beam transfer function as 



^tme 



^tme 



omeas omeas 



(39) 

where the average is over CMB fluctuations for fixed lenses. The 
statistics of ^tme^^meas discusscd in Planck Collaboration 
VII (2013), where it is found that the beam uncertainty may be 
described by a small number A^eig of eigenmodes (typically five). 



(Dmeas \ ^eig 



(40) 



Here, E^' are the uncertainty eigenmodes (indexed by /) for the 
beam determination, and gi are independent Gaussian random 
variables with unit variance. The frequency-band eigenmodes 
are small for 100-217 GHz (no more than 1% at / < 2048), and 
so we can equate 



Dtme / omeas 

^ -1-ln' ^ 



omeas 



^tme 



(41) 



We can also ignore terms of order E^ and propagate the inde- 
pendent beam uncertainties directly to a corresponding set of 
independent eigenmodes for 0. These eigenmodes are given as 



E^^ = L 



with the model that 



(42) 



(43) 



In Fig. 7 we plot the beam uncertainty eigenmodes as well 
as the corresponding normalization eigenmodes. The normal- 
ization eigenmodes are flatter than the beam uncertainty eigen- 
modes, and do not converge to zero at low-L, reflecting the fact 
that at all scales the lensing estimate takes most of its weight 
from modes of £i and £2 close to the noise/beam cutofl' (at 
^ ^ 1000 for 100 GHz and ^ ^ 1500 for 143 and 217 GHz). 
As can be seen in Fig. 7, the beam normalization eigenmode is 
essentially twice the value of the beam uncertainty eigenmode 
on these scales. 

Owing to the small size of the eigenmodes, the propagation 
to a fractional uncertainty on C^^ is accomplished essentially 
by applying a factor of two to the eigenmodes of Eq. (42). The 
normalization uncertainty due to beams for the lensing poten- 
tial power spectrum is therefore approximately 1% at 100 GHz 
and 0.5% at 143 and 217 GHz. When computing the beam un- 
certainty for our minimum- variance result, we take the 143 and 
217 GHz beam eigenmodes to be uncorrected. 

5.3. Cosmological 

Similarly to the case with the beam transfer function of the previ- 
ous subsection, uncertainty in the cosmological model may also 
be propagated to an uncertainty in the estimator normalization. 
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Fig. 7. Propagation of eigenmodes for beam transfer function 
uncertainty (left-hand panels) to an uncertainty in the lens- 
ing normalization (right-hand panels). The normalization eigen- 
modes are flat, with amplitudes given by approximately twice 
the beam uncertainty at the beam and noise cutofl' scale (at 
£ ^ 1000 for 100 GHz and £ ^ 1500 for 143 and 217 GHz). 



There are two aspects in which cosmological uncertainty can en- 
ter our analysis. First, when estimating the AC^'^lj^^ contribution 

of Eq. (15) we assume a fiducial C^^, obtained from the param- 
eters of the fiducial ACDM cosmology described in Sect. 1 . We 
take a random subset of samples from Planck MCMC chains for 
the ACDM cosmology and calculate the corresponding AC^'^lj^^ 
bias. Compared to our fiducial correction, we see small-scale 
fluctuations with an amplitude of 0.2% with a typical scale of 
AL = 150, superimposed on larger overall amplitude scatter with 
a standard deviation of approximately 6%. This result is reason- 
able: in the ACDM model the Planck TT measurements have the 
power to detect the smoothing efl'ects of gravitational lensing on 
the acoustic peaks at a significance level of lOcr, corresponding 
to a 10% constraint on the overall amplitude of C^^, consistent 
with the result from analysis of individual samples above. To ac- 
count for cosmological uncertainty in the AC^'^lj^^ correction, 
we therefore assign it a 10% overall amplitude uncertainty. This 
correction is small for the L range considered in our fiducial like- 
lihood, and so this uncertainty is essentially negligible compared 
to our statistical error bars (although we include it in our error 
budget nevertheless). 

The second point at which cosmological uncertainty enters 
our analysis is in the normalization of our estimates. In Eq. (12), 
as well as in our tests of the normalization in Sect. 4, we have 
assumed that the statistically-anisotropic covariance induced by 
lensing is equal to that assumed when weighting the estimator, 
given by Eq. (6). If the true cosmological power spectrum diff'ers 
from the fiducial one used in Eq. (6) then we will misestimate 



f>LM 



= i + aV = 



^0 ^0,tme^ 



(44) 

where wf^^^^ is given by Eq. (6), with the modification that it is 
evaluated with the "true" ensemble-average lensed CMB power 
spectrum rather than the fiducial one Cj^. We have chosen to ac- 
count for this normalization uncertainty coherently with the tem- 
perature likelihood. Given a proposed Cj^, C^j^ pair, we renor- 

maUze ^ (1 + A[^)2cf in our likelihood Eq. (23) so that 
it is directly comparable to our measurement. 



6. Results 

In this section, we present the primary results of this paper; a 
measurement of the CMB lensing potential over approximately 
70% of the sky and its associated power spectrum estimate C^^ . 
Our lensing potential map and power spectrum have several sci- 
ence implications, which we will discuss in proceeding subsec- 
tions. 

(1) The lensing power spectrum probes the matter power spec- 
trum integrated back to the last- scattering surface, and pro- 
vides an additional low redshift (z ~ 2) lever arm to 
CMB-alone parameter constraints at z = 1100. Highlights 
of the additional information that the Planck lensing like- 
lihood contributes to the parameter constraints of Planck 
Collaboration XVI (2013) are presented in Sect. 6.1. 

(2) Our lensing potential map correlates with the CMB tem- 
perature via the late-time integrated Sachs- Wolfe (ISW) ef- 
fect sourced by dark energy. We measure the cross- spectrum 
between our lens reconstruction and the low-^ temperature 
multipoles in Sect. 6.2, which we find is present at the level 
predicted in ACDM, and discrepant with the null hypothesis 
of no correlation at approximately 2. 5 a. 

(3) Our lensing potential map correlates with the fluctuations 
of the cosmic infrared background, which also trace large- 
scale structure and have good redshift overlap with the CMB 
lensing potential. We do not discuss this correlation here, but 
instead refer to the detailed analysis of Planck Collaboration 
XVIII (2013). 

(4) Finally, our lensing potential map correlates with external 
tracers of large-scale structure. In Sect. 6.3 we present cross- 
correlations with several representative galaxy catalogues. 

Our primary results are based on a minimum- variance (MV) 
combination of the 143 and 217 GHz Planck maps, although as 
an initial consistency check we will also present most of our re- 
sults in this section for the individual frequency maps as well. In 
Sect. 7 we will present a more extensive suite of consistency tests 
on our power spectrum estimates and lensing map. In particular 
in Sect. 7.1, we will compare with results obtained on difl'erent 
masks and on the foreground-cleaned Planck maps described in 
(Planck CoUaboration XII 2013). 

Our MV lensing map is shown in Fig. 8, where we plot 
the Wiener-filtered reconstruction of the lensing potential. At 
the noise levels of our reconstruction, approximately half of the 
modes in this map are noise. At a visual level, the lensing map of 
Fig. 8 has comparable signal-to-noise to the COBE-DMR maps 
of the CMB temperature anisotropy (Bennett et al. 1996). Future 
measurements, including the Planck full mission release which 
contains approximately twice the amount of data used here, will 
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improve on this first full- sky map of the CMB lensing poten- 
tial. As is illustrated in the simulated reconstruction of Fig. 4, 
there will be clear visual correlations between this map and fu- 
ture measurements. 

In Fig. 10 we plot the power spectra of our individual 100, 
143, and 217 GHz reconstructions as well as the minimum- 
variance reconstruction. The agreement of all four spectra is 
striking. Overall, our power spectrum measurement is reason- 
ably consistent with the ACDM prediction, given our measure- 
ment error bars. Dividing the L g [1, 2048] multipole range into 
bins of AL = 64 and binning uniformly in [L(L + l)]^C^^, we 
obtain a reduced x'^ for the diff'erence between our power spec- 
trum estimate and the model of 40.7 with 32 degrees of freedom. 
The associated probability to exceed is 14%. On a detailed level, 
there are some discrepancies between the shape and amplitude 
of our power spectrum and the fiducial model however. Our like- 
lihood is based on the multipole range 40 < L < 400, which 
captures 90% of the available signal-to-noise for an amplitude 
constraint on C^^. This range was chosen as the region of our 
spectrum least likely to be contaminated by systematic eff'ects 
(primarily uncertainties in the mean-field corrections at low-L, 
and uncertainties in the Gaussian and point-source bias correc- 
tions at high-L). Estimating an average amplitude for the fiducial 
lensing power spectrum for a single bin over this multipole range 
using Eq. (25) we find an amplitude of A4o^4oo = 0-94 ± 0.04 
relative to the fiducial model (which has A = 1). The power in 
this region is consistent with the fiducial model, although l.Scr 
low (the corresponding probability-to-exceed for the of this 
diff'erence is 15%). The low- and high-L extent of our likelihood 
were deliberately chosen to have enough expected lensing signal 
to enable a lOcr detection of lensing on either side, bookending 
our likelihood with two additional consistency tests. On the low- 
L side, we have a good agreement with the expected power. As 
will be discussed in Sect. 7.4, our measurement at L < 10 fails 
some consistency tests at a level comparable to the expected sig- 
nal. The L < 10 modes, which we suspect are somewhat con- 
taminated by errors in the mean-field subtraction, are neverthe- 
less consistent with the fiducial expectation, as can be seen in 
Fig. 10; we measure Ai^io = 0.44 + 0.54. Extending to the lower 
limit of our likelihood, with a single bin from 10 < L < 40 we 
measure Aio^4o = 1.02 ± 0.12. On the high-L side of our fidu- 
cial likelihood, there is tension however. Extending from the fi- 
nal likelihood multipole at L = 400 to the maximum multipole 
of our reconstruction, we find A4oo^2048 = 0-68 ±0.13, which 
is in tension with A = 1 at a level of just over 2.4cr. The rel- 
atively low power in our reconstruction is driven by a dip rel- 
ative to the ACDM model spectrum between 500 < L < 750, 
as can be seen in Fig. 10. We show this feature more clearly 
in the residual plot of Fig. 11. This deficit of power is in turn 
driven by the 143 GHz data. For an estimate of the power spec- 
trum using only 143 GHz, we measure ^Iq1^2048 = 0-37 + 0.18. 
The 217 GHz reconstruction is more consistent with the model, 
having A^q1^2048 ~ ^'^^ - These two measurements are 
in tension; we have A^Q^^^^iMs ~ ^-^^ + 0.18, which is a 2.5cr 
discrepancy. The error bar on this diff'erence accounts for the ex- 
pected correlation between the two channels due to the fact that 
they see the same CMB sky. A larger set of consistency tests 
will be presented in Sect. 7. We note for now that the bins from 
40 < L < 400 used in our likelihood pass all consistency tests, 
and show better agreement between 143 and 217 GHz. Although 
L < 40 and L > 400 are not included in our nominal likelihood, 
when discussing the use of the lensing likelihood for cosmo- 
logical parameter constraints in the following section we will 




Galactic South 



Fig. 8. Wiener-filtered lensing potential estimate 

^LM = ^L^i^LM - ^LM^ rcconstruction, in Galactic 

coordinates using orthographic projection. The reconstruction 
is bandpass filtered to L g [10,2048]. The Planck lens recon- 
struction has S/N < I for individual modes on all scales, so 
this map is noise dominated. Comparison between simulations 
of reconstructed and input in Fig. 4 show the expected level 
of visible correlation between our reconstruction and the true 
lensing potential. 
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Fig. 9. Wiener-filtered lensing potential estimates, as in Fig. 8, 
for the individual 143 and 217 GHz maps. 
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Fig. 10. Lensing potential power spectrum estimates based on the individual 100, 143, and 217 GHz sky maps, as well our fiducial 
minimum- variance (MV) reconstruction which forms the basis for the Planck lensing likelihood. The black line is for the best-fit 
ACDM model of Planck Collaboration XVI (2013). 



perform additional cross-checks using these bins to ascertain 
whether they would have any significant implications for cos- 
mology. 

In addition to the Planck power spectrum measurements, in 
Fig. 11 we have overplotted the ACT and SPT measurements 
of the lensing potential power spectrum (Das et al. 2013; van 
Engelen et al. 2012). It is clear that all are very consistent. 
The Planck measurement has the largest signal-to-noise of these 
measurements; as we have already discussed the 40 < L < 400 
lensing likelihood provides a 4% constraint on the amplitude of 
the lensing potential power spectrum, while the constraint from 
current ACT and SPT measurements are 32% and 16% respec- 
tively. These measurements are nevertheless quite complemen- 
tary. As a function of angular scale, the full-sky Planck power 
spectrum estimate has the smallest uncertainty per multipole of 
all three experiments at L < 500, at which point the additional 
small-scale modes up to ^max = 3000 used in the SPT lensing 
analysis lead to smaller error bars. The good agreement in these 
estimates of C^^ is reassuring; in addition to the fact that the ex- 
periments and analyses are completely independent, these mea- 
surements are sourced from fairly independent angular scales 
in the temperature map, with i < 1600 in the case of Planck, 
i < 2300 in the case of ACT, and ^ < 3000 in the case of SPT. 
Cross-correlation of the Planck lensing map with these indepen- 
dent measures of the lensing potential will provide an additional 
cross-check on their consistency, however at the power spectrum 
level they are already in good agreement. 

6.1. Parameters 

Weak gravitational lensing of the CMB provides sensitivity 
to cosmological parameters aff'ecting the late-time growth of 
structure which are otherwise degenerate in the primary CMB 



anisotropics imprinted around recombination. Examples include 
the dark energy density in models with spatial curvature and the 
mass of neutrinos that are light enough (my < 0.5 eV) still to 
have been relativistic at recombination. 

To connect our measurement of the lensing power spectrum 
to parameters, we construct a lensing likelihood nominally based 
on the multipole range 40 < L < 400, cut into eight equal- width 
bins with AL = 45 to maintain parameter leverage from shape 
information in addition to our overall amplitude constraint. In 
Table 1 we present bandpowers for these eight bins using the in- 
dividual 100, 143, and 217 GHz reconstructions as well as the 
MV reconstruction which is the basis for our nominal likeli- 
hood. The bandpower estimates and their uncertainties are bro- 
ken down into constituent parts as discussed in Sect. 2. Based on 
these bandpowers, we form a likelihood following Eq. (23). The 
measurement errors on each bin are measured by Monte-Carlo 
using 1000 simulations, and the bins are sufficiently wide that 
we can neglect any small covariance between them (this is dis- 
cussed further in Appendix D). We analytically marginalize over 
uncertainties that are correlated between bins, including them in 
the measurement covariance matrix. This includes beam transfer 
function uncertainties (as described in Sect 5.2), uncertainties in 
the point source correction (Sect. 7.2) and uncertainty in the N^^"^ 
correction. 

As the lensing likelihood is always used in conjunction with 
the Planck TT power spectrum likelihood, we coherently ac- 
count for uncertainty in Cj^ by renormalizing our lensing po- 
tential measurement for each sample, as described in Sect. 5.3. 

The lensing likelihood is combined with the main Planck 
TT likelihood (Planck Collaboration XV 2013) - constructed 
from the temperature (pseudo) cross-spectra between detec- 
tor sets at intermediate and high multipoles, and an exact ap- 
proach for Gaussian temperature anisotropics at low multipoles 
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100 GHz Reconstruction Bandpowers 


-^min -^max 


^RAW _ ^NO ^PS A^"*" ^MC ^BEAM 


A 


L^Cf X 10^ 


40 84 


( 7.883 - 6.948) + 0.191 +0.005 + 0.003 +0.035 + 0.003 0.007 +0.011 


0.889 + 0.191 


7.297 + 1.568 


85 129 


(11.077- 9.999) + 0.197 +0.008 + 0.004 +0.040 + 0.004 -0.010 +0.011 


1.040 + 0.197 


6.711 + 1.272 


130 174 


(16.360 - 15.411) + 0.256 +0.014 + 0.007 +0.045 + 0.004 0.025 +0.011 


0.866 + 0.257 


4.294 + 1.273 


175 219 


(25.031 - 23.896) + 0.337 +0.024 + 0.013 +0.057 + 0.006 0.043 +0.011 


1.011+0.337 


3.947 + 1.316 


220 264 


(36.642 -35.881) + 0.442 +0.040 + 0.021 +0.104 + 0.010 0.009 +0.010 


0.607 + 0.443 


1.920+ 1.400 


265 309 


(51.002 -50.170) + 0.556 +0.058 + 0.030 +0.205 + 0.020 -0.076 +0.010 


0.645 + 0.557 


1.693 + 1.462 


310 354 


(66.348 - 65.291) + 0.707 +0.071 + 0.037 +0.304 + 0.030 -0.010 +0.010 


0.692 + 0.708 


1.534 + 1.571 


355 400 


(81.994 -81.084) + 0.793 +0.077 + 0.040 +0.332 + 0.033 -0.121 +0.010 


0.621 + 0.795 


1.177 + 1.507 


40 400 


(17.725 - 16.731) + 0.112 +0.016 + 0.008 +0.062 + 0.006 0.004 +0.011 


0.912 + 0.113 








143 GHz Reconstruction Bandpowers 


-^min -^max 


^RAW _ ^NO ^PS J^l ^MC ^BEAM 


A 


L^Cf X 10^ 


40 84 


( 3.936 - 2.841) + 0.089 +0.005 + 0.002 +0.026 + 0.003 -0.021 +0.004 


1.086 + 0.090 


8.917 + 0.736 


85 129 


( 5.008 - 4.058) + 0.090 +0.007 + 0.004 +0.032 + 0.003 -0.016 +0.004 


0.928 + 0.090 


5.991 +0.582 


130 174 


( 7.076- 6.088) + 0.109 +0.012 + 0.006 +0.040 + 0.004 -0.020 +0.004 


0.956 + 0.109 


4.743 + 0.540 


175 219 


(10.092- 9.207) + 0.137 +0.021 +0.010 +0.058 + 0.006 0.013 +0.004 


0.793 + 0.138 


3.094 + 0.538 


220 264 


(14.254- 13.279) + 0.174 +0.032 + 0.016 +0.107 + 0.011 -0.010 +0.004 


0.845 + 0.175 


2.671 + 0.553 


265 309 


(19.049- 18.063) + 0.211 +0.045 + 0.022 +0.191 +0.019 0.008 +0.004 


0.742 + 0.213 


1.948 + 0.558 


310 354 


(24.524 - 22.862) + 0.244 +0.055 + 0.026 +0.271 + 0.027 -0.01 1 +0.004 


1.347 + 0.246 


2.988 + 0.547 


355 400 


(28.355 - 27.576) + 0.273 +0.060 + 0.029 +0.298 + 0.030 -0.006 +0.004 


0.427 + 0.276 


0.808 + 0.524 


40 400 


( 8.403 - 7.401) + 0.048 +0.016 + 0.008 +0.064 + 0.006 -0.013 +0.004 


0.935 + 0.050 








217 GHz Reconstruction Bandpowers 


-^min -^max 


^RAW _ ^NO ^PS J^l ^MC ^BEAM 


A 


L^Cf X 10^ 


40 84 


( 3.884 - 2.742) + 0.088 +0.003 + 0.003 +0.026 + 0.003 -0.007 +0.004 


1.120 + 0.088 


9.200 + 0.721 


85 129 


( 4.817- 3.866) + 0.088 +0.005 + 0.005 +0.032 + 0.003 -0.012 +0.004 


0.927 + 0.089 


5.982 + 0.572 


130 174 


( 6.601 - 5.788) + 0.103 +0.007 + 0.009 +0.040 + 0.004 -0.000 +0.004 


0.766 + 0.103 


3.800 + 0.512 


175 219 


( 9.544- 8.637) + 0.131 +0.012 + 0.014 +0.058 + 0.006 0.012 +0.004 


0.824 + 0.132 


3.217 + 0.516 


220 264 


(13.409- 12.317) + 0.162 +0.018 + 0.021 +0.103 + 0.010 0.004 +0.004 


0.967 + 0.163 


3.057 + 0.516 


265 309 


(17.825 - 16.583) + 0.202 +0.025 + 0.028 +0.179 + 0.018 -0.003 +0.004 


1.041 +0.204 


2.733 + 0.536 


310 354 


(22.391 - 20.915) + 0.227 +0.030 + 0.034 +0.251 + 0.025 -0.014 +0.004 


1.210 + 0.231 


2.683 + 0.513 


355 400 


(26.119 -25.068) + 0.265 +0.032 + 0.037 +0.279 + 0.028 -0.026 +0.004 


0.766 + 0.269 


1.451+0.510 


40 400 


( 8.122- 7.111) + 0.046 +0.009 + 0.011 +0.064 + 0.006 -0.004 +0.004 


0.942 + 0.048 








MV Reconstruction Bandpowers 


-^min -^max 


^RAW _ ^nO ^ps ^n1 ^MC ^BEAM 


A 


L^Cf X 10^ 


40 84 


( 3.297- 2.170) + 0.073 +0.006 + 0.003 +0.026 + 0.003 -0.011 +0.003 


1.105 + 0.073 


9.077 + 0.601 


85 129 


( 4.062 - 3.062) + 0.074 +0.009 + 0.004 +0.032 + 0.003 -0.014 +0.003 


0.972 + 0.074 


6.276 + 0.476 


130 174 


( 5.480 - 4.542) + 0.084 +0.015 + 0.007 +0.040 + 0.004 -0.012 +0.003 


0.895 + 0.085 


4.442 + 0.420 


175 219 


( 7.616- 6.758) + 0.103 +0.024 + 0.011 +0.058 + 0.006 0.009 +0.003 


0.766 + 0.104 


2.992 + 0.406 


220 264 


(10.587- 9.576) + 0.131 +0.036 + 0.017 +0.101 +0.010 -0.003 +0.003 


0.876 + 0.133 


2.771+0.419 


265 309 


(14.004- 12.903) + 0.156 +0.048 + 0.022 +0.175 + 0.017 0.000 +0.003 


0.878 + 0.158 


2.304 + 0.415 


310 354 


(17.677 - 16.285) + 0.181 +0.057 + 0.026 +0.246 + 0.025 -0.009 +0.003 


1.099 + 0.185 


2.438 + 0.409 


355 400 


(20.647 - 19.583) + 0.208 +0.061 + 0.028 +0.274 + 0.027 -0.024 +0.003 


0.753 + 0.212 


1.427 + 0.401 


40 400 


( 6.788 - 5.768) + 0.039 +0.020 + 0.009 +0.066 + 0.007 -0.009 +0.003 


0.943 + 0.040 





Table 1. Reconstruction bandpower amplitudes for 100, 143, and 217 GHz as well as the minimum- variance reconstruction, com- 
puted following Eq. (25). The overall amplitude A has been evaluated as A = A^^"^ - A^^ - A^^ - A^^ - A^^. Note that most 
of the systematic errors are significantly correlated between bands. These bandpowers are wide enough that the correlation of the 
measurement uncertainty between bins is negligible, as discussed in Appendix D. The final column lists the physical bandpower 
estimates and their errors obtained by scaling flat averages of L'^C^^ across each bin in our fiducial ACDM model by the given A 
value. 
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Fig. 11. Replotting of Fig. 10, removing 100 GHz for easier 
comparison of 143 and 217 GHz. Also plotted are the SPT band- 
powers from van Engelen et al. (2012), and the ACT bandpow- 
ers from Das et al. (2013). All three experiments are very consis- 
tent. The lower panel shows the difference between the measured 
bandpowers and the fiducial best-fit ACDM model. 



- in Planck Collaboration XVI (2013) to derive parameter con- 
straints for the six-parameter ACDM model and well-motivated 
extensions. Lensing also aff'ects the power spectrum, or 2-point 
function, of the CMB anisotropics, and this eff'ect is accounted 
for routinely in all Planck results. On the angular scales rele- 
vant for Planck, the main eff'ect is a smoothing of the acoustic 
peaks and this is detected at around lOcr in the Planck tempera- 
ture power spectrum (Planck Collaboration XVI 2013). The in- 
formation about C^^ that is contained in the lensed temperature 
power spectrum for multipoles i 1 3000 is limited to the ampli- 
tude of a single eigenmode (Smith et al. 2006). In extensions of 
ACDM with a single additional late-time parameter, lensing of 
the power spectrum itself can therefore break the geometric de- 
generacy (Stompor & Efstathiou 1999; Sherwin et al. 2011; van 
Engelen et al. 2012; Planck Collaboration XVI 2013). As dis- 
cussed in Appendix D and Schmittfull et al. (2013), cosmic vari- 
ance of the lenses produces weak correlations between the CMB 
2-point function and our estimates of C^*^, but they are small 
enough that ignoring the correlations in combining the two like- 
lihoods should produce only sub-percent underestimates of the 
errors in physical cosmological parameters. 

In the following, we illustrate the additional constraining 
power of our C^^ measurements in ACDM models and one- 
parameter extensions, highlighting those results from Planck 
Collaboration XVI (2013) where the lensing likelihood is influ- 
ential. 



6.1 .1 . Six-parameter ACDM model 

In the six-parameter ACDM model, the matter densities, Hubble 
constant and spectral index of the primordial curvature perturba- 
tions are tightly constrained by the Planck temperature power 
spectrum alone. However, in the absence of lensing the am- 
plitude As of the primordial power spectrum and the reioniza- 
tion optical depth r are degenerate, with only the combination 
As^~^^, which directly controls the amplitude of the anisotropy 
power spectrum on intermediate and small scales being well de- 
termined. This degeneracy is broken by large-angle polarization 
since the power from scattering at reionization depends on the 
combination A^r^. In this first release of Planck data, we use 
the WMAP nine-year polarization maps (Bennett et al. 2012) in 
combination with Planck temperature data. With this data com- 
bination, C^j^ is rather tightly constrained in the ACDM model 
(see Fig. 12) and the direct measurements reported here provide 
a non-trivial consistency test of the model. 

The eight C^^ bandpowers used in the lensing likelihood are 
compared to the expected spectrum in Fig. 12 (upper-left panel). 
For the latter, we have used parameter values determined from 
the main Planck likelihood in combination with WMAP polar- 
ization (hereafter denoted WP) and small-scale power spectrum 
measurements (hereafter highL) from ACT (Das et al. 2013) and 
SPT (Reichardt et al. 20 12)^ In this plot, we have renormalized 
the measurements and their error bars (rather than the theory) us- 
ing the best-fit model with a variant of the procedure described 
in Sect. 5.3. Since the lensed temperature power spectrum in the 
best-fit model is very close to that in the fiducial model used 
to normalise the power spectrum estimates throughout this pa- 
per, the power spectrum renormalisation factor (1 -h A^^)^ of 

Eq. (44) is less than 0.5% in magnitude. The predicted C^j^ in 
the best-fit model diff'ers from the fiducial model by less than 
2.5% for L < 1000. The best-fit model is a good fit to the mea- 
surements, with = 10-9 and the corresponding probability 
to exceed equal to 21%. Significantly, we see that the ACDM 
model, calibrated with the CMB fluctuations imprinted around 
z= 1 100, correctly predicts the evolution of structure and geom- 
etry at much lower redshifts. The 68% uncertainty in the ACDM 
prediction of C^^ is shown by the dashed lines in the upper-left 
panel of Fig. 12. We can assess consistency with the direct mea- 
surements, properly accounting for this uncertainty, by introduc- 
ing an additional parameter A^ that scales the theory C^^ in the 
lensing likelihood. (Note that we choose not to alter the lensing 
efl'ect in Cj^.) As reported in Planck Collaboration XVI (2013), 
we find 



Af^ = 0.99 : 



: 0.05 (68%; Planck-\-lQnsmg-\-WF-\-highL), 



in excellent agreement with Af^ = 1 . 

An alternative route to breaking the Ag-r degeneracy is pos- 



sible for the first time with Planck. Since C^^ is directly propor- 
tional to As, the lensing power spectrum measurements and the 
smoothing eff'ect of lensing in Cj^ (which at leading order varies 
as Ag ^"^^) can separately constrain Ag and r without large-angle 
polarization data. The variation of C^^ with r in ACDM models 



^ As discussed in detail in Planck Collaboration XVI (2013), the pri- 
mary role of the ACT and SPT data in these parameter fits is to constrain 
more accurately the contribution of extragalactic foregrounds which 
must be carefully modelled to interpret the Planck power spectra on 
small scales. For ACDM, the foreground parameters are sufficiently de- 
coupled from the cosmological parameters that the inclusion of the ACT 
and SPT data has very little effect on the cosmological constraints. 
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Fig. 12. Upper left: Planck measurements of the lensing power spectrum compared to the ACDM mean prediction and 68% con- 
fidence interval (dashed lines) for models fit to P/a^c^+WP+highL (see text). The eight bandpowers are those used in the Planck 
lensing likelihood; they are renormalized, along with their errors, to account for the small diff'erences between the lensed Cj^ in 
the best-fit model and the fiducial model used throughout this paper. The error bars are the ±\cr errors from the diagonal of the 
covariance matrix. The colour coding shows how C^j^ varies with the optical depth r across samples from the ACDM posterior 
distribution. Upper right: as upper-left but using only the temperature power spectrum from Planck. Lower left: as upper-left panel 
but in models with spatial curvature. The colour coding is for VIk. Lower right: as upper-left but in models with three massive 
neutrinos (of equal mass). The colour coding is for the summed neutrino mass 2 ruy. 



constrained only by the Planck temperature power spectrum is 
illustrated in the upper-right panel of Fig. 12, and suggests that 
the direct C^j^ measurements may be able to improve constraints 
on T further. This is indeed the case, as shown in Fig. 13 where 
we compare the posterior distribution of r for the Planck temper- 
ature likelihood alone with that including the lensing likelihood. 
We find 

T = 0.097 ± 0.038 (68%; Planck) 

T = 0.089 + 0.032 (68%; P/a^c^-hlensing). 

At 95% confidence, we can place a lower limit on the optical 
depth of 0.04 (Planck-\-lensmg). This very close to the optical 
depth for instantaneous reionization at z = 6, providing further 
support for reionization being an extended process. 

The T constraints via the lensing route are consistent with, 
though weaker, than those from WMAP polarization. However, 
since the latter measurement requires very aggressive cleaning 
of Galactic emission (see e.g. Fig. 17 of Page et al. 2007), the 
lensing constraints are an important cross-check. 



6.1 .2. Effect of the large and small scales on the 
six-parameter ACDM model 

Before exploring the further parameters that can be constrained 
with the lensing likelihood, we test the eff'ect on the ACDM 
model of adding the large-scale (10 < L < 40) and small-scale 
(400 < L < 2048) lensing data to our likelihood. Adding addi- 
tional data will produce random shifts in the posterior distribu- 
tions of parameters, but these should be small here since the mul- 
tipole range 40 < L < 400 is designed to capture over 90% of the 
signal-to-noise (on an amplitude measurement). If the additional 
data is expected to have little statistical power, i.e., the error bars 
on parameters do not change greatly, but its addition produces 
large shifts in the posteriors, this would be symptomatic either 
of internal tensions between the data or an incorrect model. 

In Fig. 14, we compare the posterior distributions of the 
ACDM parameters for Planck-\-WF-\-highL alone with those af- 
ter combining with various lensing likelihoods. Adding our fidu- 
cial lensing likelihood (second column) reduces the errors on pa- 
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Fig. 13. Marginalised constraints on the optical depth in ACDM 
models from the Planck temperature power spectrum (Planet, 
solid black), and additionally including the lensing likeli- 
hood (P/a^c^+lensing; dashed red) or WMAP polarization 
(Planck-\-WF', dashed-dotted blue). We use a prior r > 0.01 in 
all cases. 



rameters by a small amount and the median values shift by rather 
less than Icr for all parameters. The largest gain is for flch^ (and 
Hq) where the errors improve by 20%. Adding further large- and 
small-scale data produces no significant reduction in error bars, 
as expected. For most parameters, the medians also change very 
little except for Qch^ which is dragged low by a further O.Scr on 
adding the small-scale lensing information. (The shift in Hq is 
due to the anti-correlation between Hq and flch^ caused by the 
acoustic- scale degeneracy in the temperature power spectrum; 
see Planck Collaboration XVI 2013.) These findings are consis- 
tent with the power spectrum amplitude measurements discussed 
in Sect. 6: we can lower the lensing power by reducing the matter 
density, and this is favoured by the lower amplitudes measured 
from the small-scale lensing power spectrum. 

The tension between the small-scale power and the power 
over the L = 40-400 range included in our fiducial likelihood, 
coupled with our lower confidence in the accuracy of the bias 
removal on small scales, is the reason that we do not include 
these smaller scales at this stage in the Planck lensing likelihood. 

6.1 .3. Spatial curvature and dark energy 

Inflation models with sufficient number of ^-folds of expansion 
naturally predict that the Universe should be very close to be- 
ing spatially flat. Constraining any departures from flatness is 
therefore a critical test of inflationary cosmology. However, the 
primary CMB anisotropics alone sufl'er from a geometric degen- 
eracy, whereby models with identical primordial power spec- 
tra, physical matter densities and angular-diameter distance to 
last-scattering have almost identical power spectra (Efstathiou 
& Bond 1999). The degeneracy is partly broken by lens- 
ing (Stompor & Efstathiou 1999), with small additional con- 
tributions from the late-ISW efl'ect (on large scales) and by 
projection efl'ects in curved models (Howlett et al. 2012). In 
ACDM models with curvature, the geometric degeneracy is two- 
dimensional, involving the curvature and dark energy density, 
and this limits the precision with which either can be determined 
from the CMB alone. 
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Fig. 14. Marginalized posteriors for the six-parameter ACDM 
model, shown as box plots, for Planck-\-WF-\-highL with various 
lensing likelihoods. The red and blue lines are the median and 
mean, respectively. The box and bar correspond to 68% and 95% 
of the probability density, both centered on the median. The left- 
most column is without the lensing likelihood and the median 
of these constraints is shown by the grey line. The remaining 
columns show the efl'ect of adding in the fiducial lensing like- 
lihood (second column), and further adding a low-L bin (third 
column), high-L bins (fourth column) or both (final column). 

With the high-significance detection of lensing by Planck 
in the temperature power spectrum (Planck Collaboration XVI 
2013), and via the lens reconstruction reported here, the geomet- 
ric degeneracy is partially broken, as shown in Fig. 15. The long 
tail of closed models with low dark energy density (and expan- 
sion rate at low redshift) allowed by the geometric degeneracy 
have too much lensing power to be consistent with Planck's mea- 
sured temperature and lensing power spectra (see also Fig. 12). 
We find marginalised constraints on the curvature parameter of 



= -0.042 



+0.027 
0.018 



(68%; Planck-^WF-^highL) 



= -0.0096!^:X (68%;P/a^cy^-hlensing-hWP-hhighL), 

so that lensing reconstruction reduces the uncertainty on by 
more than a factor of two over limits driven by the smoothing 
eff'ect on the acoustic peaks of Cj^. This improvement is consis- 
tent with the spread in C^^ in curved models constrained by the 
temperature power spectrum, relative to the errors on the recon- 
struction power spectrum; see Fig. 12. Note that the mean value 
of Qk also moves towards zero with the inclusion of the C^^ 
measurements. Adding the high-L and low-L data to the likeli- 
hood brings no more than a percent-level improvement on the 
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Fig. 15. Two views of the geometric degeneracy in curved ACDM models which is partially broken by lensing. Left: the degeneracy 
in the Q^-^a plane, with samples from P/a^c^+WP+highL colour coded by the value of H^. The contours delimit the 68% and 95% 
confidence regions, showing the further improvement from including the lensing likelihood. Right: the degeneracy in the Q^k-Hq 
plane, with samples colour coded by Qa. Spatially-flat models lie along the grey dashed lines. 



constraint. We see that the CMB alone now constrains the ge- 
ometry to be flat at the percent level. Previous constraints on 
curvature via CMB lensing have been reported by SPT in com- 
bination with the WMAP-1 data: Qk = -O.OOS^^-^j^ (68%; Story 
et al. 2012). This constraint is consistent, though almost a factor 
of two weaker, than that from Planck. Tighter constraints on cur- 
vature result from combining the Planck data with other astro- 
physical data, such as baryon acoustic oscillations, as discussed 
in Planck CoUaboration XVI (2013). 

Lensing efl'ects provide evidence for dark energy from the 
CMB alone, independent of other astrophysical data (Sherwin 
et al. 201 1). In curved ACDM models, we find marginalised con- 
straints on Da of 



= 0.57 
= 0.67 



+0.073 
■0.055 
+0.027 
■0.023 



(68%; Planck+WF+highL) 

(68%; P/a^c^+lensing-FWP+highL). 



Again, lensing reconstruction improves the errors by more than 
a factor of two over those from the temperature power spectrum 
alone. 



6.1.4. Neutrino masses 

The unique eff'ect in the unlensed temperature power spectrum 
of massive neutrinos that are still relativistic at recombination 
is small. With the angular scale of the acoustic peaks fixed 
from measurements of the temperature power spectrum, neutrino 
masses increase the expansion rate at z > 1 and so suppress clus- 
tering on scales larger than the horizon size at the non-relativistic 
transition (Kaplinghat et al. 2003). This eff'ect reduces C^^ for 
L > 10 (see Fig. 12) and gives less smoothing of the acoustic 
peaks in Cj^. As discussed in Planck Collaboration XVI (2013), 
the constraint on 2 nty from the Planck temperature power spec- 
trum (and WMAP low-£ polarization) is driven by the smoothing 
eff'ect of lensing: ^m^ < 0.66 eV (95%; Planck-\-WF-\-highL). 
Curiously, this constraint is weakened by additionally including 
the lensing likelihood to 



my < 0.85 eV, (95%; P/a^cy^+WP+highL), 



reflecting mild tensions between the measured lensing and tem- 
perature power spectra, with the former preferring larger neu- 



trino masses than the latter. Possible origins of this tension are 
explored further in Planck Collaboration XVI (2013) and are 
thought to involve both the C^*^ measurements and features in 
the measured Cj^ on large scales (£ < 40) and small scales 
^ > 2000 that are not fit well by the ACDM-hforeground model. 
As regards C^^, Fisher estimates show that the bandpowers in 
the range 130 < L < 309 carry most of the statistical weight 
in determining the marginal error on 2 ^y, and Fig. 12 reveals 
a preference for high 2 from this part of the spectrum. (We 
have checked that removing the first bandpower from the lensing 
likelihood, which is the least stable to data cuts and the details 
of foreground cleaning as discussed in Sect. 7, has little impact 
on our neutrino mass constraints.) We also note that a similar 
trend for lower lensing power than the ACDM expectation on 
intermediate scales is seen in the ACT and SPT measurements 
(Fig. 11). Adding the high-L information to the likelihood weak- 
ens the constraint further, pushing the 95% limit to 1.07 eV. This 
is consistent with our small-scale measurement having a signifi- 
cantly lower amplitude. At this stage it is unclear what to make 
of this mild tension between neutrino mass constraints from the 
4-point function and those from the 2-point, and we caution 
over-interpreting the results. We expect to be able to say more 
on this issue with the further data, including polarization, that 
will be made available in future Planck data releases. 



6.2. Correlation witJn tine /SI/1/ Effect 

As CMB photons travel to us from the last scattering surface, 
the gravitational potentials that they traverse may undergo a non- 
negligible amount of evolution. This produces a net redshift or 
blueshift of the photons concerned, as they fall into and then 
escape from the evolving potentials. The overall result is a con- 
tribution to the CMB temperature anisotropy known as the late- 
time integrated Sachs-Wolfe (ISW) eff'ect, or the Rees-Sciama 
(R-S) eff'ect depending on whether the evolution of the poten- 
tials concerned is in the linear (ISW) or non-linear (R-S) regime 
of structure formation (Sachs & Wolfe 1967; Rees & Sciama 
1968). In the epoch of dark energy domination, which occurs af- 
ter z ~ 0.5 for the concordance ACDM cosmology, large-scale 
potentials tend to decay over time as space expands, resulting 
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in a net blueshifting of the CMB photons which traverse these 
potentials. 

In the concordance ACDM model, there is significant over- 
lap between the large-scale structure which sources the CMB 
lensing potential and the ISW eff'ect (greater than 90% at 
L < 100), although it should be kept in mind that we cannot 
observe the ISW component by itself, and so the eff'ective cor- 
relation with the total CMB temperature is much smaller, on the 
order of 20%. 

The correlation between the lensing potential and the ISW 
eff'ect results in a non-zero bispectrum or three-point function 
for the observed CMB fluctuations. This bispectrum is peaked 
for "squeezed" configurations, in which one short leg at low-^ 
supported by the ISW contribution is matched to the lensing- 
induced correlation between two small-scale modes at high- 
i. Constraints on the amplitude of the lensing-ISW bispec- 
trum using several diff'erent estimators are presented in Planck 
Collaboration XXIV (2013). Here we will present an additional 
constraint, in which the bispectrum measurement is recast as an 
estimate for the amplitude of the cross-spectrum C^*^, using the 
filtering and frequency map combinations of our baseline lensing 
reconstruction. Our measurements are in good agreement with 
those made in Planck Collaboration XXIV (2013); a detailed 
comparison of several lensing-ISW bispectrum estimators, in- 
cluding the one used here, is presented in Planck Collaboration 
XIX (2013). 

Following Lewis et al. (201 1), we begin with an estimator for 
the cross-spectrum of the lensing potential and the ISW eff'ect as 



C 



T(p _ -^sky , 

~ + 1 V 

M 



(45) 



where T^m = Q T^m is the Wiener-filtered temperature map 
and ^ is given in Eq. (13). In Fig. 16 we plot the measured cross- 
spectra for our individual frequency reconstructions at 100, 143, 
and 217 GHz as well as the MV reconstruction. We also plot the 
mean and scatter expected in the fiducial ACDM model. 

To compare quantitatively the overall level of the measured 
C^^ correlation to the value in ACDM, we estimate an overall 
amplitude for the cross-spectrum as 



= f^T^ £ (2L + l)C[^-'^'^-CJ^(C[^Arf ). (46) 

Li—J^mm 

The overall normalization N^^ is determined from Monte-Carlo 
simulations. For our processing of the data, we find that it is well 
approximated (at the 5% level) by the analytical approximation 



£ (2L+1)(C[^'^^-) /(C[^A^f) 



(47) 



The estimator above is equivalent to the KSW and skew-Q es- 
timators of Komatsu et al. (2005); Munshi et al. (2011b) for the 
lensing-ISW bispectrum which are used in Planck Collaboration 
XXIV (2013) (up to implementation details such as filtering). 
The mean-field subtraction performed when computing 0lm can 
be identified with the linear term of Creminelli et al. (2006) 
which is necessary to minimize the estimator variance. The con- 
tribution to the total S /N of this estimator as a function of the 
short leg L is plotted in Fig. 2, where it can be seen that the con- 
straining power for the fiducial correlation is almost entirely at 
L < 100. 
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Fig. 16. Lensing-ISW bispectrum-related cross spectra com- 
puted from Eq. (45). Black dashed lines indicate the average 
value for simulations, while dark/light gray filled regions indi- 
cate the expected one/two standard deviation scatter, also mea- 
sured from simulations. The thin magenta line gives the expected 
C^*^ cross- spectrum for our fiducial model. The agreement of 
this curve with the simulation average illustrates that our esti- 
mator is accurately normalized. In all the quantitative analysis 
of this section we ignore L < 10, although we have plotted the 
cross- spectra at these multipoles for interest. 



In Table 2 we present measured values for the amplitude of 
the lensing-ISW bispectrum using Eq. (45). The uncertainties on 
A^*^ are determined by Monte-Carlo. We use the multipole range 
10 < L < 100, given some of the potential systematic issues 
with these multipoles identified in Sect. 7.4, although as can be 
seen from Fig. 16, the inclusion of lower multipoles does not sig- 
nificantly aff'ect our results. Note that for the ISW-lensing mea- 
surements, inaccuracies in the mean-field subtraction do not bias 
the estimator although they may degrade the statistical errors on 
large scales. The differences between the diff'erent amplitude fits 
are well within the expected scatter, as we show in Table 3. 

As a point of interest, we have also split our amplitude 
constraint into the contribution from even and odd multipoles. 
There are well known odd/even-multipole power asymmetries in 
the temperature anisotropics on large angular scales, the study 
of which is somewhat limited by the small number of avail- 
able modes (Land & Magueijo 2005; Kim & Naselsky 2010; 
Gruppuso et al. 2011; Bennett et al. 2011). The lensing potential 
gives a potentially new window on these power asymmetries, as 
a third somewhat independent measurement of power on large 
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Table 2. Best-fit amplitudes of the lensing-ISW bispectrum A^^ 
given by Eq. (46) with Lmin = 10 and L^ax = 100, as well as 
split into even/odd-L contributions. All of the fits are consistent 
(within 2cr) with the amplitude of unity expected in our fiducial 
ACDM model. 



Lensing-ISW Amplitudes 




A^^ (all L) 


A^"^ (even L) 


A^^ (odd L) 


100 GHz . . 
143 GHz . . 
217 GHz . . 
MV 


. . 0.93 + 0.52 
. . 0.81+0.36 
. . 0.87 + 0.35 
. . 0.78 + 0.32 


0.45 + 0.72 
0.27 + 0.48 
0.54 + 0.49 
0.25 + 0.45 


1.44 + 0.73 
1.37 + 0.52 
1.22 + 0.49 
1.32 + 0.46 



Table 3. Diff'erences between the lensing-ISW amplitude fits 
of Table 2, along with Icr scatter determined from Monte-Carlo 
simulations. 





Lensing 


^-ISW Pair Differences 




143 GHz 


217 GHz 


MV 


100 GHz . . . 


. +0.13 + 0.38 


+0.06 + 0.42 


+0.16 + 0.41 


143 GHz . . . 




-0.07 + 0.26 


+0.03 + 0.18 


217 GHz . . . 






+0.10 + 0.13 


Even- Odd . . 


. -1.10 + 0.69 


-0.68 + 0.69 


-1.07 + 0.64 



angular scales. As we can see in Table 2, there is a large diff'er- 
ence between the odd/even-L contributions in our lensing-ISW 
bispectrum estimate, related to the odd/even temperature asym- 
metry. 

In Table 3 we present diff'erences between the amplitude 
measurements of Table 2, as well as the expected scatter, ac- 
counting for correlations between the estimates due to common 
CMB and (in the case of MV vs. 143 or 217 GHz) noise fluc- 
tuations. We see that our estimates are all very consistent. The 
diff'erence in measured odd/even multipole power, while strik- 
ing, is not statistically significant at greater than 2cr. 

To conclude, we see a correlation between our lensing poten- 
tial measurement and the large-scale temperature anisotropics 
that is consistent with level expected due to the present-day 
dark energy domination in our fiducial ACDM cosmology. Our 
amplitude measurements are in agreement with those presented 
in Planck Collaboration XXIV (2013); these independent mea- 
surements are compared in detail in Planck Collaboration XIX 
(2013). 

6.3. Correlation with Galaxy Catalogues 

The CMB lensing potential is an integrated measure of all mat- 
ter in the Universe back to the last scattering surface, with a 
broad kernel which peaks at z ~ 2, but has significant contri- 
butions from both lower and higher redshifts. There are there- 
fore significant correlations between the CMB lensing potential 
and other tracers of large-scale structure. Such correlations have 
already been observed between lensing potential from WMAP 
and the NVSS (Smith et al. 2007; Hirata et al. 2008, on the order 
of 3 (T in significance), as well as between the lensing potential 
from SPT with optical and infrared catalogs from BCS, WISE, 
and Spitzer/IRAC (Bleem et al. 2012, between 4 and 5cr), and 
the ACT lensing potential with quasars from SDSS (Sherwin 
et al. 2012, at 3.8(t). Here we will show several representative 
examples of correlations with the Planck lensing potential, be- 



tween 7cr and 20cr in significance. Our goal here is not to per- 
form an in depth study, but rather to demonstrate the power 
of our public lensing map. In addition to the correlations pre- 
sented here, there is a powerful correlation between the Cosmic 
Infrared Background (CIB). Correlation of the Plancklcnsing 
with the CIB fluctuations as probed by the highest frequency 
Planck channels is observed at greater than 40cr in significance, 
and has been subjected to a more detailed analysis and mod- 
elling, which is presented in an accompanying paper, Planck 
CoUaboration XVIII (2013). 

To predict expected levels of correlation for a given galaxy 
catalog, we use the Limber approximation (Limber 1954) with a 
simple linear bias model, in which the cross- spectrum between 
two mass tracers is given by 



J dxKl(x)4(x)P(k = LIx.xX (48) 



where ;^ is conformal distance and P(k,x) is the 3D matter power 
spectrum for wavenumber k at conformal lookback time x- 
and are kernels associated with lensing and with the galaxy 
catalog of interest respectively, and are given by 



L2 



(49) 



Here b(z) is a redshift-dependent linear bias parameter, and 
dN/dz describes the redshift distribution of the galaxy popula- 
tion. As with the lensing-ISW bispectrum, we can construct a 
simple pseudo-C^ estimator for the cross correlation using 



2L 



skyj 1* 



(50) 



where /skyj is the sky fraction common to both the lens recon- 
struction and the catalog, and giM is the harmonic transform of 
the galaxy fractional overdensity. Denoting the positions of the 
objects p in the catalog as hp, the transform is given by 



(51) 



where N^^ is the surface density of objects in gal/steradian. 
With a fiducial model C^*^'^^' for cross spectrum obtained using 
Eq. (48), the minimum-variance estimator for its overall ampli- 
tude is 



(2L+l)Cf^'-Cf 



4'(Cf + A^^^)(Cf + 



L 



(52) 



where C^^ is the signal power spectrum of the catalog, which 
can be estimated using Eq. (48) with for both weight func- 
tions. Here we have defined the spectrum A^j^ as a scaling at 
each multipole L of the cross-correlation power spectrum. The 
normalization W^*^ is given by 



2 



(2L+l)(cf'^'^-f 



(53) 



In Fig. 17 we plot the contributions to A^'^ as a function of 
L for several surveys which have significant correlations with 
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Fig. 17. Cross-spectra of the Planck MV lensing potential with several galaxy catalogs, scaled by the signal-to-noise weighting 
factor A^*^ defined in Eq. (52). Cross-correlations are detected at approximately 20cr significance for NVSS, lOcr for SDSS LRGs 
and la for both MaxBCG and WISE. 



the Planck MV lensing potential: the NVSS quasar catalog, the 
MaxBCG cluster catalog, an SDSS LRG catalog, and an infrared 
catalog from the WISE satellite. The error bars for each correla- 
tion are measured from the scatter of simulated lens reconstruc- 
tions correlated with each catalog map, and are in generally good 
agreement (at the 20% level) with analytical expectations. These 
catalogs are discussed in more detail below. 

1. NVSS Quasars: The NRAO VLA Sky Survey (NVSS) 
(Condon et al. 1998) is a catalog of approximately two mil- 
lion sources north of ^ = -40° which is 50% complete at 
2.5mJy. Most of the bright sources are AGN-powered ra- 
dio galaxies and quasars. We process this catalog follow- 
ing Smith et al. (2007), pixelizing the catalog at HEALPix 
A^side = 256 and projecting out the azimuthally symmetric 
modes of the galaxy distribution in ecliptic coordinates to 
avoid systematic striping eff'ects in the NVSS dataset. We 
model the expected cross-correlation for this catalog using 
a constant b{z) =1.1 and a redshift distribution centered at 
zo = 1.1 given by 



dN 



exp(- 
exp(- 



{z-zof \ 
0.8)2 ) 



2(0.^. 
{z-zof \ 
2(0.3)2 ) 



{z < Zo) 

iz > zo). 



(54) 



For this model, in the correlation with the MV lens recon- 
struction we measure an amplitude of A^y^^ = 1.03 ± 0.05. 
2. SDSS LRGs: We use the LRG catalog of Ross et al. (201 1); 
Ho et al. (2012) based on Sloan Digital Sky Survey Data 
Release 8 (SDSS DR8), which covers 25% of the sky. After 
cutting to select all sources with photometric redshift 0.4 < 
z < 0.8, and /7gai > 0.2, we are left with approximately 1.4 x 
10^ objects with a mean redshift of z = 0.55 and a scatter 
of ±0.07. Apart from the cut above, we do not perform any 
additional weighting on /7gai. We model this catalog using 
dN/dz taken from the histogram of photometric redshifts. 



and take b(z) = 2. We measure A^^^ = 0.96 ± 0.10, very 
consistent with expectation. 

3. MaxBCG Clusters: The MaxBCG cluster catalog (Koester 
et al. 2007) is a collection of 13, 823 clusters over approx- 
imately 20% of the sky selected from the SDSS photomet- 
ric data, covering a redshift range 0.1 < z < 0.3. It is be- 
lieved to be 90% pure and more than 85% complete for 
clusters with M > 1 x lO^^M©. To simplify the sky cover- 
age, we have discarded the three southern SDSS stripes in 
the catalog, which reduces the overall sky coverage to ap- 
proximately 17%. There are accurate photometric redshifts 
(A^ ~ 0.01) for all objects in the catalog, and so we can 
construct dN/dz directly from the histogram of the redshift 
distribution. Although these clusters are at very low red- 
shift compared to the typical structures which source the 
CMB lensing potential, they are strong tracers of dark mat- 
ter, with an eff'ective bias parameter of b(z) = 3 (Huetsi 
2009). We obtain a similar average bias parameter {b(M, z)} 
for the MaxBCG clusters if we combine the mass-richness 
relation of Bauer et al. (2012) and the halo bias prescription 
of Tinker et al. (2010). Here measure a correlation with the 
Planck lensing potential of ^MaxBCG ~ 1.54 ± 0.21. This 
is significantly larger than expected given the simple model 
above, although as can be seen in Fig. 17 the shape of the 
correlation is reasonable agreement. 

4. WISE Catalog: The Wide Field Survey Infrared Explorer 
(WISE) satellite (Wright et al. 2010) has mapped the full 
sky in four frequency bands Wl — W4 at 3.4, 4.6, 12, and 
22yum respectively. We start from the full mission catalog, 
which contains over five hundred and sixty million objects. 
To obtain a catalog with roughly uniform sensitivity over 
the full sky and to eliminate stellar contamination we fol- 
low Kovacs et al. (2013), selecting all sources with Wl mag- 
nitudes less than 15.2 at galactic latitudes greater than 10°, 
and require Wl - W2 > 0.2 and W2 - W3 > 2.9. We cut 
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all sources which are flagged as potentially spurious, or for 
which more than 30% of frames observing a given source 
are marked as contaminated by moon-glow. These cuts leave 
us with approximately 2,308,751 sources. We additionally 
remove all sources which lie outside an /sky = 0.6 Galaxy 
mask. Following Goto et al. (2012), we model the redshift 
distribution of these sources using 



dN 



oc zexp 



_2 \ 



(55) 



with zo = 0.1, cr^ = 0.1, and where the proportionality con- 
stant is chosen such that J dzdN/dz = 1 . The bias is taken 
to be Z? = 1. For this model we measure an amplitude of 
"^wisE ~ ^-^^ ±0.13. Note that the cuts we have used above 
are fairly conservative- in particular, relaxing the magnitude 
cut above we can obtain much larger cross-correlation sig- 
nals (on the order of 40cr in significance). 



It can be seen from Fig. 17 that we have significant correla- 
tions with the objects in all four of these catalogues, with ampli- 
tudes and shapes roughly at the expected level. The correlations 
here are only a demonstrative sample of what is possible with 
the Planck lensing potential map. Note that although the lens- 
ing potential has a peak sensitivity to high-redshift tracers (such 
as NVSS), significant correlations may also be observed with 
relatively low redshift catalogs. From the perspective of cross- 
correlation, the Planck lensing potential represents a mass sur- 
vey with a well understood redshift kernel over the full-sky, and 
we anticipate that correlations with it will prove useful as a direct 
probe of the correlation between luminous and dark matter. 



7, Consistency and Systematic Tests 

In this section we perform a number of consistency and system- 
atic tests. Most of our tests will focus on the robustness of the 
MV reconstruction which forms the basis for our fiducial lensing 
likelihood. In Fig. 18 we summarize several of our main consis- 
tency checks: foreground tests obtained by comparing individ- 
ual frequencies, adjusting the Galactic mask level, or comparing 
to results obtained analyzing component separated maps (dis- 
cussed in Sect. 7.1), tests of our correction at the power spec- 
trum level for contamination by unresolved point- sources (in 
Sect. 7.2), tests for possible noise bias obtained by taking ap- 
propriate correlations of maps with independent noise realiza- 
tions (discussed further in Sect. 7.3), and tests using alternative 
estimators which are bias hardened against various systematics 
(Sect. 7.4). In all cases, we find that our MV power spectrum is 
generally consistent with these alternative reconstructions. 

Our baseline reconstruction method uses the methodology 
outlined in Sect. 2, as well as a specific choice of filtering to 
produce the T which are the inputs to our lensing estimators. We 
test our robustness to these choices in Sect. 7.5, as well as our 
implementation of our main reconstruction pipeline, by compar- 
ing to three alternative and largely independent pipelines. 

In addition to these targeted tests, we have also performed 
several generic null tests for unexpected signals in our data. 
These are deferred to Appendix A, where we study for survey- 
to- survey discrepancies, lensing curl modes, and the map-level 
statistics of our lens reconstruction. 



7.1. Foreground Contamination 

In this section, we test the consistency of our results with respect 
to the treatment of Galactic foreground contamination. We use 
several diff'erent tests to gain perspective on possible issues: 

(1) We compare individual frequency reconstructions. The 
agreement of the lensing power spectra estimated at 100, 
143, and 217 GHz is apparent within the error bars plot- 
ted in Fig. 10. However, because a large fraction of the un- 
certainty in these estimates comes from CMB fluctuations 
themselves which are common between frequencies they are 
significantly correlated. In Fig. 18 we plot the diff'erence be- 
tween the 143 and 217 GHz reconstructions (the scatter for 
100 GHz is too large to provide a useful test of the MV re- 
construction so we have not plotted it, although it is also 
consistent). Both are quite consistent, although we can see 
that the somewhat large amplitude for our first likelihood 
bin, just above L = 40 is being driven by 217 GHz. 

(2) We perform lens reconstruction using both more aggressive 
and more conservative Galaxy masks than our fiducial /sky = 
0.7 analysis mask, constructed following the description in 
Sec. 3. The results of this test are also shown in Fig. 18, for 
masks with /sky = 0.6 and /ky = 0.8. We can see a highly 
significant discrepancy at very low-L between our fiducial 
/]^y = 0.7 mask and the more aggressive /ky = 0.8 mask, 
however other bins (in particular in the 40 < L < 400 range 
are untouched. Our fiducial results are completely consistent 
with the more conservative /ky = 0.6 results. 

(3) We perform lens reconstruction on the component-separated 
CMB maps from Planck Collaboration XII (2013). These 
are the product of four diff'erent foreground-removal algo- 
rithms: RULER, NiLC, SEVEM, and SMicA. Each of these methods 
combines the full set of nine Planck frequency maps from 
30 857 GHz to obtain a best-estimate CMB map. This 
makes a detailed characterization of these maps more diffi- 
cult. However analysis of these maps does provide a valu- 
able consistency test. In Fig. 18 we compare the reconstruc- 
tions obtained on these component separated maps to our 
MV reconstruction, analyzed using the same Galaxy mask. 
The foreground-cleaned maps are generally in good agree- 
ment with our fiducial reconstruction, although the first bin 
above L = 40 is consistently about one standard-deviation 
lower than our MV result, which as we have already pointed 
out is driven high by 217 GHz. The confidence masks pro- 
duced in Planck Collaboration XII (2013) cover significantly 
more sky than the /ky = 0.7 mask used in our baseline anal- 
ysis. 

The remainder of this subsection describes our analysis of the 
component separated maps in more detail. For the power spec- 
trum tests of Fig. 18, we analyze these maps following the same 
methodology as for our fiducial results. As a mask, we use the 
union of our fiducial /ky = 0.7 Galaxy mask, point source masks 
at 100, 143, and 217 GHz, as well as the "method" masks pro- 
vided for each component separation algorithm. We character- 
ize the bias terms of Eq. (15) as well as the statistical uncer- 
tainties of the lens reconstruction on each map using the same 
lensed CMB and FFP6 noise realizations which are discussed 
in Sect. 4. These have been run through each of the compo- 
nent separation algorithms above, resulting in a set of simulated 
foreground-cleaned maps (note that no foreground signal was 
included in the simulations which were run through the sepa- 
ration algorithms, and so these simulations only reproduce the 
mixing employed by each algorithm of the individual frequency 
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Fig. 18. Summary of internal consistency tests between our fiducial minimum- variance (MV) reconstruction and a set of alternatives 
designed to test sensitivity to potential issues. The top panel shows C^^ estimates, with measurement error bars. The bottom panels 
show the residual with respect to the MV reconstruction in units of the MV measurement uncertainty. The gray band marks the Icr 
deviation uncertainty of the MV reconstruction. The error bar on each data point in the lower panels gives the standard deviation 
of the scatter between each result and the MV, determined from Monte Carlo simulations which account for the correlated CMB, 
noise and foreground power between estimators. Comparison of the uncertainty on the scatter points and the gray band gives an 
indication of how constraining each test is. The various tests are described in more detail in subsections of Sect. 7. 
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maps). To match the power spectrum of these simulations to the 
power spectrum of the data maps, we find it is necessary to add 
extragalactic foreground power following the model in Sect. 4, 
with Acib = 18 fiYs? and Asrc = 28 yu^^. The resulting simula- 
tions have a power spectrum which agrees with that of the CMB 
map estimate based on the data to better than 2% at / < 2048. 
This could be improved slightly by tailoring a specific correc- 
tion for each map. We also add homogeneous pixel noise with a 
level of 12 yuK arcmin. If we neglected this power, the agreement 
would be only at the 8% level, primarily due to the noise term 
(the Acib and Asrc contributions are each at the level of 1 - 2%). 
Due to the procedure which we use to subtract the disconnected 
noise bias (Eq. 17) from our lensing power spectrum estimates, 
the inclusion of these components does not significantly aff'ect 
our results, but comparison with the values used for our single- 
frequency simulations in Sect. 4 are a useful indicator of the ex- 
tent to which the foreground separation algorithms are able to 
remove extragalactic foreground power in the high-^ regime. 

As already discussed, our results on the component- 
separated CMB maps are presented in Fig. 18. Because the 
CMB and FFP6 noise components of the foreground-cleaned 
map simulations are the same as those used to characterize 
our fiducial lens reconstruction, we can measure the expected 
scatter between the foreground separated maps and our fidu- 
cial reconstruction. This scatter will be slightly overestimated 
because we have not attempted to coherently model the con- 
tribution to the reconstruction noise from residual diff'use ex- 
tragalactic foreground power. For the eight bins in 40 < L < 
400 on which our fiducial likelihood is based, we measure a 

for the diff'erence between our fiducial reconstruction and 
the corresponding foreground-cleaned reconstruction of - 
(3.14,4.3,2.5,14.7) for nilc, smica, sevem, and ruler respec- 
tively. These values associated have probability-to-exceed 
(PTE) values of (79%, 64%, 86%, 2%) respectively. At the level 
which we are able to test, the nilc, smica, and sevem foreground- 
cleaned maps give results which are quantitatively consistent 
with our fiducial reconstruction. There is more scatter between 
our fiducial reconstruction and the ruler map than expected 
from simulations, as evidenced by a very high x^ for the dif- 
ference, however as can be seen in Fig. 18, there are not any 
clear systematic diff'erences. Indeed, the discrepancy for the bins 
plotted in Fig. 18 (which diff'er somewhat from the linear bins 
used in our likelihood) is much less significant than for the bins 
of our fiducial likelihood. 

When using the component separated maps above, we have 
used the same /sky = 0.7 Galactic mask as for our MV result, al- 
though the confidence regions associated with each foreground 
cleaned map allow more sky, ranging up to /sky = 0.94 for the 
NILC method. We have used the metis pipeline (described later 
in Sect. 7.5) to test whether this improved sky coverage could 
benefit our lens reconstruction. The same method has been used 
in (Planck Collaboration XII 2013) to evaluate possible biases 
to lens reconstruction induced by these methods using the FFP6 
simulated CMB realization, described in Planck Collaboration I 
(2013), indicating that the diff'erent component separation algo- 
rithms do not alter significantly the lensing signal (at the level 
which can be tested on a single simulation). Analyzing the nilc 
map, which has the largest confidence region, we find that we 
can increase the usable sky surface up to /ky = 0.87 without 
encountering significant Galactic contamination. In Fig. 19 we 
show the striking improvement in sky coverage on the nilc map. 
SMicA and sevem are very similar; we have not considered ruler 
because of its larger noise level. 



Power spectrum estimates at this mask level show consis- 
tency with the MV reconstruction within two standard devia- 
tions of the measurement uncertainty. The increased sky cover- 
age does not bring significant improvements in the error-bars of 
the power spectrum, however. Using Eq. 20 as an estimate of the 
power spectrum variance, the larger sky coverage yields only a 
3.5% improvement at L < 40 over the MV result, decreasing 
down to at L = 400. This could be due to the diff'erent weight- 
ing used in the component separation compared to the one of 
the MV map, which results in slightly noisier maps for our pur- 
pose. While the component separated maps allow for a reduced 
mask maintaining a robust lensing potential estimation, they lead 
to a marginal improvement of the power spectrum uncertainties. 
Nevertheless, their agreement with the MV result is reassuring. 




Fig. 19. Wiener-filtered potential maps in Galactic coordinates, 
as in Fig. 8, plotted here in Mollweide projection. Top is the MV 
reconstruction, bottom is an extended reconstruction on the nilc 
component- separated map. 



7.2. Point Source Correction 

As can be seen in Table 1, the unresolved point source shot 
noise correction in any individual band for our MV likelihood 
is on the order of a few percent, reaching up to 6% for the 
highest multipole bands. Averaged over the 40 < L < 400 
band, the shot noise correction amounts to a 2% shift in the am- 
plitude of C^^, which is small but non-negligible compared to 
our statistical uncertainty of 4%. Physically, the amplitude of 
our source corrections are reasonable; at 143 GHz we measure 
^143 ~ (1.3 ± 0.6) X 10~^^ /uK^. From the radio point source 
model of De Zotti et al. (2010), this corresponds to an eff'ec- 
tive flux cut of approximately 150mJy at this frequency, roughly 
comparable to that expected for the S/N > 5 cut we make when 
masking sources in our fiducial analysis (Planck Collaboration 
XXVIII 2013). The shot noise measured at 217 GHz is lower, as 
expected given the smaller contribution from radio sources, with 
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S\^^ = (0.4 ± 0.4) X lO-^V^"^- The shot noise level which we 
measure for the MV reconstruction lies between these two, with 

= (0.6 ± 0.3) X 10-1 V^"^. 
As a consistency test of our S/N>5 source masks, we con- 
struct two new sets of source masks which only exclude S /N > 
10 sources (to check the stability of our reconstruction under 
significantly higher source contamination), or which exclude all 
detected sources with S/N> 4.2 (to test the conservatism of our 
fiducial flux cut). The results of this test are shown in the masks 
panel of Fig. 18. Although not shown explicitly in the figure, the 
source correction term with the S /N > 10 mask is considerably 
larger than our fiducial correction. For this mask we measure 

= (2.8 ± 0.3) X 10-1 V^"^. now a 9cr detection of shot noise 
in the trispectrum. The importance of the correction grows cor- 
respondingly, becoming a 9% overall correction to the C^*^ am- 
plitude in 40 < L < 400. With the source correction, we measure 
^^^400 = O-^^ ± 0-04, consistent with the A^^4qq = 0.94 ± 0.04 
measurement for our fiducial source mask. Without the source 
correction we would have measured A^Y,4oo ~ 1-04 ± 0.04, sig- 
nificantly discrepant. We can see that even for this extremely 
permissive source mask, our correction for the point source shot- 
noise trispectrum does a good job of rendering C^^ insensitive to 
source contamination. Turning to the S /N > 4.2 mask, we mea- 
sure ^MV,SNR4 = (-^-1 ± 0-3) X 10-lV^^ reduced from a 2cr 
measurement of shot noise to one no very consistent with zero. 
For this reconstruction we measure A4o^4oo = 0-93 ± 0.04, in 
good agreement with the MV result, although there does ap- 
pear to be an overall trend for lower power above L = 400 
in the S /N > 4.2 mask. At a detailed level, both of these tests 
do fail x'^ t^st on the diff'erence with the MV reconstruction. 
As the change in sky fraction between these three reconstruc- 
tions is small (the diff'erence between the fiducial 3.nd S /N>4.2 
mask is 1%, and even smaller for the S /N > 10 mask), even the 
small diff'erences quoted above are significantly larger than they 
should be if the source correction were perfect. This can be seen 
in Fig. 18, where the error bars on the diff'erence between each 
of these reconstructions and the MV is too small to be visible 
for the plotted points. For the 40 < L < 400 range on which our 
likelihood is based, we believe that the discrepancy between the 
power spectra estimated with these source masks is acceptable. 
Particularly given the results of the 5'/A^ > 10 cut, the source cor- 
rection which we perform appears to behave reasonably, and re- 
duces the sensitivity of our estimates to unmasked point sources. 
In the following section, we will also use a modified lensing esti- 
mator which has zero response to the shot noise trispectrum and 
therefore does not require a source correction. As with the tests 
above, we will find good overall agreement with the standard 
estimator used for the MV reconstruction. 

7.3. Instrumental Noise Bias 

Errors in our instrumental noise model are a source of concern 
for our lensing power spectrum estimates, in which all four legs 
of the trispectrum estimate come from a single map. In this sit- 
uation there are three Wick contractions of the noise which con- 
tribute to the power spectrum estimate in two distinct ways, de- 
pending on whether they couple between the two quadratic esti- 
mators in the power spectrum estimate or not: 

(1) The Wick contraction which couples noise contributions in- 
side the same estimator leads to a mean-field term as in 
Eq. (9). This term is sourced by the statistical anisotropy 
of the noise (due to the uneven hit distribution across the 
sky, or the correlations along the scan direction induced by 



deconvolution of the bolometer time constant). The shape 
of this mean-field is dictated by the Planck scan strategy, 
which is coherent over large scales, and so this mean-field is 
primarily at low-L, below the Lmin = 40 cutoff' of our likeli- 
hood. It can be avoided by forming a quadratic estimator in 
which the two input maps have diff'erent noise realizations. 
We perform a null test against this contribution in Fig. 18, 
where the (143x2 17) x( 143x2 17) panel shows a comparison 
of our MV power spectrum estimate to that obtained using 
a quadratic estimator in which one leg comes from 143 GHz 
and one from 217 GHz. We see generally good agreement 
with the MV result. For the 8 bins in 40 < L < 400 of our 
fiducial likelihood, we find slx^ for the diff'erence (account- 
ing for correlations between the two estimates) of 4.2, which 
has a PTE of 83%. The noise mean-field contribution can 
also be reduced by using an appropriate bias-hardened esti- 
mator. We will perform additional tests with such estimators 
in Sect. 7.4. 

(2) The two Wick contractions which couple noise terms across- 
estimators contribute to the disconnected bias AC^^^I^^^. This 
term describes the noise contribution to the reconstruction 
variance. Even given small infidelities in our noise simula- 
tions used to compute the disconnected noise bias, our use of 
the data covariance matrix in Eq. (17). It can be avoided by 
taking the cross-spectrum of two quadratic estimators with 
independent noise. We show the result of such a correlation 
compared to our MV reconstruction in Fig. 18, where the 
(143 X 143) X (143 x 143) panel shows the cross- spectrum 
between the individual 143 and 217 GHz reconstructions. 
Again, we see generally good agreement with the MV result. 
For the 8 bins in 40 < L < 400 of our fiducial likelihood, 
we find a for the diff'erence (accounting for correlations 
between the two estimates) of 6.5, which has a PTE of 59%. 

Based on these tests, we believe that noise bias cannot be a sig- 
nificant contaminant for our power spectrum estimates. 

7.4. Bias Hardened Estimators 

Most of the results in this paper use the standard weight func- 
tion W^^^^^ to form our quadratic estimators for the lensing po- 
tential. This weight function is a matched filter for the lensing 
eff'ects, and results in lensing estimates with minimal variance. 
As we discuss in Appendix C, however, the standard estima- 
tor has large mean-field contributions at low-L, primary due to 
masking and noise inhomogeneity. As pointed out in Namikawa 
et al. (2012a), these mean-fields can be mitigated with the use 
of "bias-hardened" lensing estimators, which use weight func- 
tions W^^^^^ specially constructed to project out these eff'ects. 
Following the notation of Appendix C, we consider an uncer- 
tain bias field zlm with associated weight function W^^^^^ which 
gives a contribution to the off'-diagonal elements of the CMB 
covariance matrix (and therefore the estimator mean-field) of 

A(r,„,r,,„,> = (^'^ i\ ^ ) wi,^,z,M. (56) 

LM ^ ' 

We can then form an estimator which is insensitive to zlm simply 
as 

The quadratic estimator (0-z) which uses the weight function 
above has a significantly reduced mean-field contribution from 
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Zlm- This bias-hardening procedure may be repeated iteratively 
to produce weight functions which are insensitive to several 
sources of mean-field simultaneously. In Fig. 20 we plot a com- 
parison at low-L between estimates with the standard estimator 
and estimators which are bias hardened against noise and mask 
mean-fields. We see generally good agreement between these 
sets of estimates for L > 10. Below L = 10, the pseudo-spectrum 
of the diff'erence between the standard and bias-hardened esti- 
mators is unexpectedly large. For this reason, we band- limit our 
tensing potential map to L > 10, and do not consider L < 10 for 
any of the quantitative results in this paper. 

In addition to low-L consistency tests, we may also use these 
bias-hardened estimators to estimate C^^ . The resulting power 
spectra are plotted in Fig. 18. Taking the diff'erence with the 8 
bins in 40 < L < 400 of our fiducial MV reconstruction, with 
scatter estimated from simulations, we find;^^ values of 6.6 and 
4.3 for the mask/noise-hardened estimator respectively. These 
have corresponding PTE values of 58% and 83% The error bars 
on C^^ obtained with the bias hardened estimators used here are 
generally between ten and twenty percent larger than the error 
bars of the standard estimator. We also construct an estimator 
which is bias hardened against the point source weight function 
^A^2^ of Eq. (2.4). This estimator has the distinction of having 
zero response to point source shot noise. Again, we find consis- 
tent results with the standard estimator, with = ^ ^ PTE 
of 76%. 



7.5. Alternative Methods 

All of the primary results in this paper use a lens reconstruc- 
tion pipeline based on the methodology outlined in Sect. 2. As 
a robustness test, both of this methodology and of its implemen- 
tation, we have implemented three completely independent re- 
construction pipelines. These independent pipelines make sig- 
nificantly diff'erent choices than our baseline approach, primar- 
ily in the choice of the inverse- variance filter function and the 
calculation of the correction terms in Eq. (15). Our main result 
of this section is Fig. 21, where we compare the lens reconstruc- 
tion power spectra for these alternative pipelines to our baseline 
results at 143 and 217 GHz, as well as the MV reconstruction. 
The agreement is excellent. The only aspect of implementation 
common to these four results (apart from the data which was 
analyzed) is the set of simulations, described in Sect. 4 used to 
estimate the various bias terms and characterize the estimator 
scatter. 

We now proceed to describe these alternative pipelines in 
somewhat more detail. We note first two common aspects of all 
three alternative pipelines, which diff'ers from our baseline re- 
sults: rather than accounting for the Galaxy and point-source 
mask in the filter function as is done for our baseline method, 
these methods use a completely diagonal filter, which is the same 
as the form (given in Eq. 10) which our baseline filtering asymp- 
totes to far from the mask boundaries. Instead, as we will dis- 
cuss, the alternative pipelines deal with the mask using diff'er- 
ent choices of data preparation and selection. A second diff'er- 
ence common to all three pipelines with our baseline approach 
is that when computing the disconnected noise bias AC^'^Ij^q, 
these methods do not use the two-point expression of Eq. (17), 
but rather the approximation to it based on Eq. (21), evaluated 
using an estimate of the data power spectrum. Because this cor- 
rection is quite large, the agreement of results which use this 
alternative method to calculate it is reassuring. We now proceed 
to describe the individual methods in more detail; the common 
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Fig. 20. Bias-hardened estimator consistency tests, for our MV 
reconstruction. Following the discussion in Sect. 7.4 we form 
estimators which are bias-hardened against the mean-fields due 
inhomogeneous noise levels ((p-N, cyan), masking (0-M, ma- 
genta), and both of these eff'ects simultaneously ((p- M - N, 
)range). The top panel shows the power spectra of the mean-field 
for each of these estimators. The middle panel shows the power 
spectra of the reconstructions themselves, compared to the aver- 
age expected from signal -hnoise (dashed black). The lower panel 
shows the power spectrum of the diff'erence between each bias- 
hardened estimator and the standard result. Dashed black lines 
give the expected average of this diff'erence (measured on simu- 
lations). We can see that the diff'erences are consistent with the 
expected scatter for all estimators for L > 10. 
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Fig. 21. Comparison of alternative lensing pipelines. The base- 
line results use the methodology of Sect. 2. Boxes are for the 
MV reconstruction, circles show the 143 GHz results, triangles, 
the 217 GHz ones. The two bottom panels show differences rel- 
ative to the MV result for 143 and 217 GHz. 

motivation in the development of each has been the reduction 
of the sharp gradients induced when masking, which can induce 
a mean-field several orders of magnitude larger than the lens- 
ing signal at low multipoles (as discussed in Appendix C). Each 
method takes a diff'erent approach to mitigating this mask eff'ect, 
as discussed below: 

(1) The method iso consists of applying the standard quadratic 
lensing estimator to the sky map after multiplying by an 
apodized Galactic mask, and filling point source holes us- 
ing local constrained Gaussian realizations of the CMB sig- 
nal -Fnoise. The mask mean-field is proportional to the power 
spectrum of the mask, and so as apodization smooths the 
mask boundary (suppressing its power spectrum on small 
scales), it correspondingly reduces the mean-field signifi- 
cantly. The combination of apodization and source filling 
makes this estimator very fast to apply to simulations but 
does require an involved set of correction terms and /sky fac- 
tors. Our implementation and calculation of this method is 
described in detail in Benoit-Levy et al. (2013). 

(2) The METIS method consists of inpainting the Galactic mask 
as well as the point source holes, using the sparse-inpainting 
algorithm described in (Abrial et al. (2007, 2008)). The re- 
sulting map resembles a full- sky CMB map, and therefore 
has no mask mean-field contribution. The inhomogeneous 



noise and beam-induced mean-fields do still have to be cor- 
rected however. Our implementation is based on that de- 
scribed in Perotto et al. (2010), with several improvements. 
In Perotto et al. (2010) lens reconstruction was performed on 
the inpainted map and then analyzed on the full-sky, how- 
ever further inspection has revealed that there are some spu- 
rious features in the lens reconstruction, localized to the in- 
painted region inside the Galactic mask. This is likely due to 
the inhomogeneous noise in Planck that was ignored in pre- 
vious work and cannot be well reproduced by the inpainter. 
We therefore remask the full-sky lens reconstruction with 
an apodized Galactic mask (as in Eq. 14) to remove these 
regions from our analysis. We follow the same procedure 
when evaluating the analytical expression for the AC^'^Ij^q 
bias, prewhitening and then applying an apodized Galactic 
mask to the inpainted temperature multipoles to estimate 
their power spectrum. Small residual biases are corrected us- 
ing the same AC^'^lj^^ procedure used in the main method. 
(3) The PATCHES method avoids the Galactic mask completely 
by cutting the sky into a collection of 410 small overlap- 
ping 10° X 10° patches centered on the locations of A^side = 8 
HEALPix pixels, which are then analyzed under the flat- sky 
approximation. Our implementation of this method is de- 
scribed in Plaszczynski et al. (2012). As with the iso method, 
point source holes are filled using constrained Gaussian re- 
alizations. The patches are extracted from a pre- whitened 
CMB map, and apodized with a Kaiser-Bessel window func- 
tion. The Fourier modes in each patch are fitted in real space 
using a fast Fourier-Toeplitz algorithm. No mean-field cor- 
rection is applied. Residual biases due to noise inhomogene- 
ity are removed using a AC^'^lj^^ correction which is found 
to be small. The patches method has been particularly useful 
in the early stages of our analysis, to identify outliers caused 
by unmasked point sources. 

As can be seen in Fig. 21, all three of these methods are in good 
agreement with the results of our baseline method, providing re- 
assurance that our results are insensitive to the precise details of 
our data filtering and reconstruction methodology. 



8. Conclusions 

The Planck maps have unprecedented sensitivity to gravitational 
lensing eff'ects. We see significant and consistent measurements 
of lensing for each of the high-resolution CMB channels at 100, 
143 and 217 GHz. Even the noisiest channel which we have con- 
sidered, 100 GHz, provides a lOcr detection of lensing, which is 
greater than all previous detections. Our fiducial lens reconstruc- 
tion, based on a minimum- variance combination of the 143 and 
217 GHz channels does even better, with a detection of lensing 
(relative to the null hypothesis of no lensing) at a significance 
of greater than 25(T. Notably, the noise on our reconstruction is 
low enough that it is no longer the limiting source of noise for 
many correlations with large-scale structure catalogs (several ex- 
amples of which we have given in Sect. 6.3). This marks a shift 
for CMB lensing, from the detection regime into that of standard 
cosmological probe. Our lensing potential map is publicly avail- 
able, and we look forward to the uses which may be found for 
it. 

The percent-level Planck lensing potential measurement 
pushes into the realm of precision cosmology, and requires care- 
ful validation tests which we have performed in Sect. 7. Our 
fiducial likelihood, based on the 40 < L < 400 range which is 
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most sensitive to lensing, passes all of the tests that we have per- 
formed at an acceptable level. Most importantly, it is consistent 
with individual frequency reconstructions, with more aggressive 
or more conservative masking, and with results obtained on more 
rigorously component separated maps. Our measurement of the 
lensing potential power spectrum C^j^ is broadly consistent with 
the ACDM expectations of Planck Collaboration XVI (2013), al- 
though there are some shape and amplitude tensions which lead 
to smaller than forecasted improvements on parameters such as 
the sum of neutrino masses My. As we saw in Sect. 6.1, our 
estimate of the potential power spectrum has direct cosmolog- 
ical implications, where it significantly improves CMB power- 
spectrum only constraints on curvature (or alternatively, dark 
energy). Our lensing likelihood also breaks the As - r degener- 
acy, allowing a measurement of the optical depth to reionization 
which is a useful independent cross-check on the values deter- 
mined with large angular scale polarization measurements. 

CMB lensing science has undergone a remarkable develop- 
ment in the last two years (Das et al. 2011; Sherwin et al. 2011; 
van Engelen et al. 2012). The next year is expected to bring even 
greater improvements. The 2500deg^ lensing survey of the South 
Pole Telescope is expected have comparable statistical power to 
the one reported here, with a highly complementary weighting in 
L- space which is more sensitive to smaller angular scale lenses. 
There is also additional Planck data, which will improve on the 
measurement reported here considerably. The results of this pa- 
per are based on the nominal mission Planck maps, collected 
during the first fifteen months of science operations. The cryo- 
gens used to maintain the 100 mK cooling stage for the Planck 
HFI instrument lasted considerably longer than the nominal mis- 
sion lifetime, to a full mission duration of approximately thirty 
months. This additional data, set to be released in 2014, results 
in map noise levels which are a factor two lower in power than 
those used here. As instrumental noise constitutes approximately 
half of our error budget, we expect a corresponding increase of 
approximately 25% in the variance of our lensing potential map 
as well as the uncertainties on C^*^. The Planck polarization data, 
also set for release in 2014, will provide an additional powerful 
probe of lensing. 
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Appendix A: Null Tests 

In this appendix we present several generic null tests of our lens 
reconstruction. In Appendix A. 1 we perform tests using lens re- 
construction on the two individual 6-month full- sky surveys con- 
tained in the first Planck data release. In Appendix. A.2 we per- 
form the curl-mode null test which is standard practice in Galaxy 
lensing. Finally, in Appendix. A. 3 we will perform several tests 
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Fig. A.l. Survey consistency null tests at 143 GHz. Solid black 
points show diff'erence of C^j^ estimated separately on the two 
six month surveys of the nominal mission. Magenta points show 
the auto-spectrum of obtained from maps made by taking the 
half-diff'erence of the two surveys (SI - S2) x (SI - S2). Yellow 
points show the cross- spectrum of obtained on half-diff'erence 
and half-sum maps (SI - S2) x (SI -h S2). The scatter on the tests 
which contain diff'erence maps is very small because (as can be 
seen in Fig. 6) most of the lens reconstruction "noise" is due to 
primary CMB fluctuations, which is nulled in these tests. 



on the isotropy and Gaussianity of our MV lensing potential 
map. 

A.1. Survey Consistency 

The Planck nominal mission consists of approximately two full- 
sky 6-month surveys. Our fiducial analysis is based on maps 
made from the full data range, however we may also analyze 
them first two surveys separately to exaggerate our sensitivity 
to eff'ects such as beam asymmetry, correlated noise along the 
scan direction, and zodiacal light. In Fig. A.l we plot several 
consistency tests exploiting the breakdown of the data into sep- 
arate surveys. We test for gross contamination of either survey 
by estimating the lensing power spectrum on each individually 
and then diff'erence the resulting C^, finding results consistent 
with the scatter expected from simple signal -hnoise simulations. 
To test the noise model more thoroughly, we also run lens recon- 
struction on a map constructed by taking the half-diff'erence of 
the two surveys, filtering it in the same way as the full nominal 
mission map. As can be seen in Fig. A.l, the power spectrum 
of this reconstruction null test is tiny compared to the measure- 
ment errors on C^^. Finally, to isolate the eff'ect of systematics 
which couple to the CMB and which flip signs between surveys 
(such as odd moments of the beam asymmetry), we take the 
cross- spectrum of lens reconstruction estimates obtained from 
the half-difl'erence of the two survey maps with those obtained 
from the half-sum. Again, this null test shows residual which is 
negligibly small compared to our measurement errors. 

A.2. Curl Mode Reconstruction 

In addition to the standard gradient lensing potential cp sourced 
by large-scale structure, we may also use the Planck maps to 
reconstruct a curl mode lensing field i//. The weight function as- 
sociated with the curl field is given by (Namikawa et al. 2012b) 

= - ^^^^^^^^^ V^ZTTkkaTT) 

xQ'.'( '"^"f^'^" )('l' -i) + (A^^2). (A.1) 
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Fig.A.2. Lensing curl power spectrum estimates for the MV 
reconstruction, as well as the individual 143 and 217 GHz fre- 
quency reconstructions. 



Using this weight function, we may construct an estimator for 
the curl-mode power spectrum, completely analogous to that for 
0. The response function of the curl estimator to point short shot- 
noise is zero, and so we neglect this correction. Consistency of 
the reconstructed curl mode with zero is a standard systematic 
test in weak lensing, as the curl field due to gravitational lens- 
ing should be negligible (Hirata & Seljak 2003b). One caveat 
for this test in the CMB lensing context is that the variance of 
the curl mode reconstruction is still aff'ected by the presence of 
gravitational lenses, by means of an 'W^^^"-type (Kesden et al. 
2003) contraction in the trispectrum (van Engelen et al. 2012). 
Our curl-mode power spectrum estimates are plotted in Fig. A. 2, 
after correction for the "M^^"-type bias using Eq. (18). This bias 
is also given explicitly for the curl mode in Benoit-Levy et al. 
(2013). We do not see evidence for anomalous power in the curl 
reconstruction, although we do find in Sect. C.4 that the curl es- 
timate is contaminated by pixelization eff'ects at a level too small 
to be seen in this plot. 

A.3. Statistical tests at the map level 

In this section we present several sanity tests of our lensing esti- 
mates at the map level. 

1 -point PDF of 0. In the left panel of Fig. A.3 we plot the his- 
togram, or 1 -point function, of our unnormalized potential es- 
timate for the MV reconstruction, bandlimited to 10 < L < 
2048. As can be seen from the histogram, the map is slightly 
non-Gaussian, as might be expected for an estimator which is 
based on the product of two Gaussian random fields. In the right 
panel of Fig. A.3 we plot the same histogram for a simulated 
MV reconstruction, which shows a comparable level of non- 
Gaussian structure. In Appendix D we show that this small non- 
Gaussianity is negligible for our lensing likelihood. 

Power Spectra on Patches To test for gross contamination of 
our reconstruction from any localized region of the sky, we use 
the PATCHES (described in Sect 7.5) to build an estimate of the 
lensing power spectra on a collection of approximately four hun- 
dred 10° X 10° tiles. We fill point source holes with local con- 
strained Gaussian realizations of the CMB signal. The patches 
were chosen to have no overlap with the Galaxy mask. 

The distribution of the power spectra estimated for each 
patch can be used as a diagnostic for deficiency in the mask and 
dust correction. Significantly contaminated but localized regions 
would appear as outliers. This procedure successfully detected 
several point sources which were overlooked in early versions 
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Fig. A.3. Histogram of pixel values for the unnormalized lensing 
estimate 0(w) for the MV reconstruction, in units of the standard 
deviation of the map. A Gaussian with cr = 1 is overplotted in 
black. Note that the y-axis is logarithmic. The left column shows 
the histogram for the data itself, while the right column shows a 
simulated reconstruction. The reconstructed map is slightly non- 
Gaussian, but in a way which is expected from simulations given 
the non-Gaussian nature of our two-point estimator. 




Fig. A.4. Bi-dimensional deflection power spectrum obtained by 
stacking the individual Fourier amplitude obtained on patches, 
after AC^'^^ subtraction and smoothing for the 143 and 
217 GHz channels. 
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Fig. A.5. Density plot of the individual deflection power spectra 
constructed on approximately four hundred 10° x 10° patches for 
the 143 GHz and 217 GHz channels. The lines are the theoretical 
1 and 2cr error for this patch size. No strong outliers are seen. 



30 



Planck Collaboration: Gravitational lensing by large-scale structures with Planck 



of our source masks. In Fig. A.5 we present results of this anal- 
ysis using our final Galaxy and point source masks for 143 GHz 
and 217 GHz. 



Isotropy on Patches. We can also use the same data to build 
a coarse test of the isotropy of the reconstructed lensing map. 
We perform lens reconstruction on patches which are oriented 
along the local spherical basis (eg, ^0), and then stack the corre- 
sponding 2D power spectrum estimates to build a hi- dimensional 
power spectrum for the 143 GHz and 217 GHz maps. This is 
shown in Fig. A.4, where we have additionally applied smooth- 
ing based on an undecimated wavelet transform to reduce the 
visual noise. The overplotted contour lines are placed at levels 
given by the azimuthally averaged power spectrum for L= 100, 
200, 300 and 400. 



Appendix B: Filtering 

In this section, we discuss our filtering methodology, which 
takes us from a set of maps if^^'^ indexed by frequency v and 
pixel p to the multipoles T^m which are fed into the quadratic es- 
timators used to measure lensing and possible contaminants. For 
an observed CMB which consists of Gaussian signal and noise, 
the optimal filtering (in the minimum-variance sense) for the 
construction of any quadratic anisotropy estimator is so-called 
"C"^" filtering, in which the observed map is beam-deconvolved 
and then divided by a signal-Fnoise covariance matrix. For the 
purposes of lens reconstruction, the Gaussian assumption above 
is quite accurate for Planck. The Planck map noise is well ap- 
proximated as Gaussian (Planck Collaboration VI 2013). The 
observed CMB signal is less-obviously well approximated as 
Gaussian, with the non-Gaussianity due to lensing detected at 
high significance, however for temperature lens reconstruction 
at Planck noise levels the diff'erence in signal-to-noise using a 
filter constructed in the Gaussian limit versus a more compli- 
cated "delensing" filter which properly accounts for lensing non- 
Gaussianity has been shown to be negligible (Hirata & Seljak 
2003a). Other sources of non-Gaussianity, such as unresolved 
point sources, are completely subdominant to either primary 
CMB or noise at all of the multipoles probed by our lens re- 
construction, and so we neglect them in the construction of our 
filters. 

For a single map, with a symmetric beam transfer function, 
full-sky coverage and homogeneous map noise levels, the 
filter is given by 
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where Y^mif^p) gives the value of the spherical harmonic at the 
center of pixel p, Bi gives the beam transfer function and Fi is 
given by 

Fi = — ^, (B.2) 
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where cr is the map noise level in units of yuK arcmin. The maps 
we use must always be masked to remove bright Galactic fore- 
grounds and point sources, and so we do not ever use the full- sky 
filter above, although as we will discuss the filter which we do 
use asymptotes to this form in regions far from the mask, and 
the simple form Fi of the filter function given in Eq. (B.2) will 
prove useful in many analytical calculations for our analysis. 



For the more realistic case of inhomogeneous map noise lev- 
els, with pixel- space noise correlations and a sky cut, the con- 
struction of is more involved. To account for the sky cut, we 
can take the noise covariance to be infinite for masked pixels. 
Reformulating the C"^ filter using the inverse of the noise co- 
variance matrix A^^^, in pixel space (defined to be zero if either 
p or p' is a masked pixel), and allowing for the joint-filtering of 
multiple input maps (with no noise correlations between them) 
we have 
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where v indexes the maps which are being combined and the 
pointing matrix J/^^ = B^^Y^mi^p) gives the value of the spher- 
ical harmonic at the center of pixel p, convolved with the ap- 
propriate beam-hpixel transfer function 5^. The matrix Cem,e'm' is 
given by 
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The filter of Eq. (B.3) is the one which we will use, with a 
slightly suboptimal choice for the noise covariance matrix Npp^ 
which is described below. 

The Planck map noise has several sources of pixel-pixel cor- 
relations for p p'.On small scales, there are correlations along 
the scan direction induced by the deconvolution of the bolometer 
time response (Planck Collaboration VII 2013). On large scales, 
there are correlations between scan rings, due to residual l/f 
noise fluctuations not cancelled by the baseline subtraction pro- 
cedure (Planck Collaboration VI 2013). Accounting for these 
pixel correlations would be computationally intractable, as the 
full covariance matrix Npp^ (an N^^^ sized object) is too large 
to work with on a modern computer, and so we neglect them 
in our noise filter. The remaining structure in the noise matrix 
is due to the variable coverage in diff'erent regions of the sky. 
Noise levels are lowest near the Ecliptic poles, which are visited 
more frequently than the region around the Ecliptic plane due 
to the Planck scan strategy. Although the variable noise level is 
tractable to include in our filtering procedure (and is indeed in- 
cluded in the C"^ filtering of e.g. Planck Collaboration XXIV 
(2013)), we have chosen to neglect this aspect of the noise, us- 
ing instead a fixed eff'ective noise level. We use an inverse-noise 
map in pixel space 



-i,v _ ( 180x60xM^ 

TT X (Tv 



4n 



pix 



(B.5) 



where Mp is a sky mask and CTy is an average map noise level 
in units of yuK arcmin. In addition to this noise weighting, we 
also project out the five modes corresponding to the monopole, 
dipole, and 857GHz Planck map (as a simple dust template). The 
noise level for each of these modes is taken to infinity using the 
Woodbury formula; putting each of these modes into a map Tp 
indexed by /, our total noise covariance matrix is given in terms 
of the inverse noise map above by 



^pp'Np 



l,v 



Z 



(B.6) 



where Tu is the overlap matrix between templates given by 



(B.7) 
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Our motivation for taking a fixed noise level is that with this ap- 
proach, in regions sufficiently far from the mask boundary, our 
filter asymptotes to the diagonal form of Eq. (B.l). This means 
that the normalization of our lensing estimates can be well- 
approximated analytically, which is very useful for the propa- 
gation of systematic eff'ects, and also that the normalization of 
our lensing estimates does not vary across the sky with noise 
level, which simplifies cross-correlation analysis. Our fil- 
ter is therefore optimally accounting for masking eff'ects, but not 
for noise correlations and inhomogeneity. We estimate the sub- 
optimality of neglecting these noise properties by calculating the 
quantity 



{SINf 
(S/NT 



^OPT^OPT 
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where F^^^ is the optimal filter and F^'^ is the suboptimal filter 
which we have actually used. This equation gives the S IN loss 
as a function of lens multipole L, however in practice we find 
that the L dependence is small enough that it suffices to quote a 
single average loss. To estimate the degradation due to ignoring 
noise correlations we set 



77OPT 



1 



(B.9) 



where Nj^^ is the power spectrum of the map noise. We find 
that the degradation due to neglect of noise correlations is small; 
less than 2% for all L < 2048 at 100 GHz, and less than 0.1% 
at 143 and 217 GHz. To calculate the degradation due to ignor- 
ing noise inhomogeneity, we determine the map noise level in 
the 3072 regions corresponding to A^side = 16 HEALPix pix- 
els, take F^^^ using Eq. (B.2) with the local noise level, and 
estimate a resulting S /N degradation using Eq. (B.8). The ne- 
glect of noise inhomogeneity is the dominant suboptimality of 
our filtering, although it is still small. We find an average S /N 
loss (averaged over the entire sky) of approximately 4% at 100, 
143, and 217 GHz, consistent with the simulation-based results 
of Hanson et al. (2009). We take this loss as justified, given the 
simpler normalization properties of our lensing estimates when 
neglecting variations in the map noise level. 



In our analysis, we account for most of these eff'ects with a cor- 
rective mean-field term, given by Eq. (9), which is determined 
using Monte Carlo simulations. In this appendix, we will break 
this mean-field down into its constituent parts and discuss each 
in more detail. As an overview of the results in this section, 
in Fig. C.l we plot estimate for the three largest mean-fields, 
due to masking, noise inhomogeneity, and beam asymmetry at 
143 GHz (100 and 217 GHz are qualitatively similar). These 
mean-fields all have most of their contributions on very large 
scales, dictated by the coherency of the scan strategy in the case 
of beam asymmetry and noise inhomogeneity, and of the large- 
scale nature of the Galactic foregrounds in the case of the sky 
mask. 




Fig. C.l. Analytical estimates for the power spectra of the largest 
low-L mean-fields 143 GHz. The various components are dis- 
cussed in more detail in Sect. C.l (mask). Sect. C.2 (noise), and 
Sect. C.3 (beams). The mean-fields all couple most strongly to 
even modes of the lens reconstruction, due to the approximate 
north/south symmetry of the scan strategy and Galactic mask. 



Appendix C: Mean-Fields 

As discussed in Sect. 2, the quadratic lensing estimators which 
we use are designed to detect statistical anisotropy induced by 
lensing. There are a number of non-lensing sources of statistical 
anisotropy which can mimic the lensing signal to some extent. 
In our analysis, the eff'ects we consider are 

(1) The application of a sky mask, which introduces sharp gra- 
dients that may be misinterpreted as lensing. 

(2) Noise inhomogeneity, which causes the overall power to 
fluctuate across the sky and can resemble the convergence 
component of lensing. 

(3) Beam asymmetry, which smears the fluctuations more along 
one direction than another and can mimic the shear compo- 
nent of lensing. 

(4) Pixelization, in which detector samples are accumulated into 
pixels, introduces a spurious deflection field on the pixel 
scale because the centroid of the hit distribution in each pixel 
does not necessarily lie at the pixel center. 



Our discussion will focus on constructing simple models for 
each source of mean-field. Following Hanson et al. (2010), we 
will identify each of the individual contributions to mean-field 
with a tracer zlm that sources a contribution to the CMB covari- 
ance matrix given by 



LM ^ 



^2 ^ \ wz 

m2 M I ^^hhL^ 



(C.l) 



where V^^^^^l ^ weight function describing how zlm couples 
multipoles. Such a contaminant leads to a bias for the standard 
lensing estimator given by 



IMF 



LM. 



(C.2) 



where the response function "R^^ is defined in Eq. (12). The ana- 
lytical forms for the mean-fields which we present here are used 
in Sect. 7.4 to construct "bias hardened" estimators which have 
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smaller mean-field corrections, as proposed by Namikawa et al. 
(2012a). 

C.1. Mask 

The application of a mask M(h) to the beam-convolved sky, fol- 
lowed by subsequent beam deconvolution, which can be mod- 
elled as a statistical anisotropy "M", with Mlm = / dhYlj^M(h) 
and weight function (Namikawa et al. 2012a) 
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The mask mean-field using our fiducial Galaxy and point source 
analysis mask (described in Sect. 3) is plotted in Fig. C.l. At 
very low-L, the mean field induced by this mask can be much 
larger than the reconstruction noise level. At L = 2 for example, 
it can be up to 10"^ times greater than the reconstruction noise 
level, although the amplitude of the mean-field is a very steep 
function of L. 

The mask mean-field is proportional to C^^, and so with 
percent-level agreement between the power spectrum of our sim- 
ulations we use to determine the mean-field and the data power 
spectrum (discussed in Sect. 4), we can expect errors on the 
power spectra of the mask mean-field at the level of one part 
in 10"^. At L = 2 for 143 and 217 GHz this corresponds to a 
level of [L(L -h ly^C"^^ ^ 4 x 10"^ comparable to the lensing 
potential power spectrum. There are several good ways to mit- 
igate the eff'ect of the mask mean-field. By apodizing the mask 
before lens reconstruction, the gradients at the mask boundary 
can be greatly reduced which significantly suppresses the mean- 
field (van Engelen et al. 2012; Benoit-Levy et al. 2013). At 
the map level, the only issue with such an apodization for our 
purposes is that it results in a position-dependent normalization 
which can complicate cross -correlation analyses. A related ap- 
proach, which we have taken in our power spectrum analysis, 
is to apodize the lens reconstruction after performing the re- 
construction. This amounts to removing the regions of near 
the mask boundary where the mean-field has all of its features. 
This "avoidance" approach can also be used easily in cross- 
correlation analyses, as in Hirata et al. (2004). A third approach 
is that of Namikawa et al. (2012a), who shows that it is possi- 
ble to construct "bias -hardened" lensing estimators which have 
smaller mean-field contributions. We have already showed con- 
sistency tests based on this approach in Sect. 7.4. In the low- 
L bias-hardened comparisons of Fig. 20, we saw discrepancies 
with the standard estimator of [L(L ly^C"^^ 4 x 10"^ and 
7 X 10"^, precisely the level expected above. 

For our publicly released lensing map, we have chosen sim- 
ply to high-pass filter our map to L > 10, above which we be- 
lieve that mean-field errors are subdominant to the reconstruc- 
tion noise. 

C.2. Instrumental Noise 

The instrumental noise of Planck is rich in statistically 
anisotropic structure, sourced by (Planck Collaboration VI 2013) 

(1) Striping eff'ects, due to residual noise correlations on large 
time scales after destriping. 



(2) Correlations along the scan direction induced by the de- 
convolution of the detector time response and the low-pass 
rolloff' filter. 

(3) Large variations in the pixel noise level due to the larger 
density of pixel hits near the Ecliptic poles, due to the scan 
strategy. 

All of these eff'ects can be avoided by forming estimators us- 
ing pairs of data with independent noise realizations, as we have 
done in Sect. 7.3, although our baseline results do not use such 
cross-correlations. Here we discuss the size of the noise mean- 
field. Items (1) and (2) are difficult to study analytically. For 
white noise with a variance in each pixel given by N(hp), item 
(3) can be identified with a statistical anisotropy "A/^" of defined 
by (Hanson et al. 2009) 
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Note that the shape of this bias is dictated by the hit count dis- 
tribution, which is known quite well, rather than the noise level, 
which is more uncertain. The shape of the noise mean-field is 
plotted in Fig. C.l. As with the mask, it is large compared to 
our signal and noise but falls off quickly with multipole. For 
L < 10, given that there are known issues with our noise mod- 
elling at the 5% level in power (as discussed in Sect. 4), it can 
be important construct quadratic estimators from correlations of 
maps with independent noise realizations (Apart from our bias- 
hardened estimator tests in Sect. 7.4, this is also done for L = 1 
in (Planck Collaboration XXVII 2013)), however at the L > 10 
multipoles of our lensing map and power spectrum, uncertainties 
in the noise mean-field should be negligible. 

C.3. Beam Asymmetry 

Planck observes the sky after convolution with a "scanning 
beam", which captures its eflfective response to the sky as a 
function of displacement from the nominal pointing direction. 
Asymmetry in the shape of this beam (for example ellipticity) 
means that the ffuctuations in the map are smeared depending 
on the beam orientation, which can mimic the shearing effects 
of gravitational lensing. Here we present semi- analytical esti- 
mates of this beam asymmetry mean-field, and discuss how sim- 
ulations of beam asymmetry are included in our simulations to 
correct for it. As with the noise inhomogeneity, we will see that 
the large-scale coherency of the Planck scan strategy means that 
beam asymmetries only mimic large-scale modes of (piM- 

Decomposing the scanning beam into harmonic coefficients 
B^fn, each time-ordered data (TOD) sample can be modelled as 
(neglecting the contribution from instrumental noise, which is 
independent of beam asymmetry) 



(C.5) 



where the TOD samples are indexed by /, and T^m is the under- 
lying sky signal. The spin spherical harmonic ^7^^^ rotates the 
scanning beam to the pointing location {6, 0), while the e~^^^' 
factor gives it the correct orientation. On the small scales which 
lensing is sensitive to, it is a good approximation to assume that 
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the procedure of mapmaking from TOD samples is just a process 
of binning: 

T{hp) = ^ TJHp, (C.6) 

iep 

where hp is the centroid of pixel p, and the sum is taken over all 
hits assigned to pixel p. Hp is the number of hits in pixel p. Here 
we will assume that the pointing for each observation of pixel 
p is identical to the pixel center h. The purpose of this assump- 
tion is simply to isolate the effect of beam asymmetry from the 
(small) effect of pixelization. We will discuss the pixelization 
effect in Sect. C.4. Placing all hits at pixel centers, combining 
Eq. (C.5) and Eq. (C.6) we have 



T(h) = ^ w(h, -s) 
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where the "scan strategy maps" w(w, s) are given by 



w(hp, -^) = ^ ( 
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where h is the center of pixel p and in the second line we have 
made a gradient approximation to Yimin) about the pixel center, 
and used the effective deflection vector 
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where hi is the pointing for observation /, which is assigned to 
pixel p (with center at hp and number of hits Hp). Implicit in this 
model is an assumption that on the arcminute pixelization scale 
the beam-convolved sky fluctuations are quite smooth, and can 
be modelled as a gradient. Pixelization therefore introduces an 
effective deflection field, given by d^^^, which biases our recon- 
struction of the lensing potential. We could in principle correct 
for the spurious lensing field which this introduces by includ- 
ing the remapping of Eq. (C.l 1) in the simulations used to esti- 
mate the mean-field correction, however as we will see it is quite 
small. We can decompose the spin-1 deflection field as 



lEp 



J d^hjl^ih) [df + id^""] {h\ (C.13) 



The scan strategy maps are spin-s objects, with corresponding 
harmonic decomposition 

s^LM = J d^hsYl^{h)w{h, s). (C.9) 

We will use Eq. (C.7) to include the effect of asymmetric beams 
in the simulations which we use to debias our lensing estimates. 
We can also propagate the effect of beam asymmetry to the lens- 
ing mean-field directly. Each mode s of the scan strategy can be 
identified with a statistical anisotropy of the type defined in 
Eq. (C.l), with = ,wlm and 
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In Fig. C.l we have plotted the expected full-sky bias for the 
best-fit Planck beams, which are measured from planet observa- 
tions in Planck Collaboration VII (2013). There are considerable 
differences in the mean-field between frequencies, as well as the 
steep scale dependence, which is dictated by the large-scale co- 
herency of the scan strategy. 

C.4. Pixelization 

In the process of mapmaking, TOD samples are accumulated 
into bins defined by HEALPix pixels. The Planck HFI maps are 
generated at HEALPix resolution 11, with a typical size of 1.'7, 
which is sufficiently close to the ^ 2.4' scale of deflections ex- 
pected from gravitational lensing effects that this remapping of 
the observed CMB is a source of concern. To model this effect, 
we begin from Eq. (C.5), but treat the beam as symmetric to 
isolate the effect of pixelization from that of beam asymmetry. 
Binning hits into pixels, we then have 
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where de and d^ are the components of d^^^ along the Q and di- 
rections respectively. The associated weight function, following 
Eq. (C.l), is given by 
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4n 
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The pixelization deflection field can also be usefully split into 
gradient- and curl-type parts, as 



(C.15) 



where the gradient and curl parts both satisfy the reality condi- 
tion 

Only the gradient mode couples to the lensing estimator, while 
the curl part only couples to the curl-mode lensing estimator 
(discussed in Sect. A.2). The statistical anisotropy introduced by 
jG,pix almost completely degenerate with that of CMB lensing, 
differing only by the fact that it operates on the beam-convolved 
CMB sky rather than the primordial one. In Fig. C.2 we plot the 
expected bias to the lensing power spectrum estimator, which is 
given by 



AC? 
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as well as the analogous bias AC^*^! . for the curl null-test es- 
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timator. To test for the presence of the pixelization lensing field 
in our reconstructions, we cross-correlate the deflection fields 

0. In Fig. C.2 we plot 



jG/c, pix ^.^j^ lensing estimates 
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and similarly for 0. We detect both d^^ and at greater than 



andJf^ 

15cr significance for all of our three main frequency bands, at the 
expected level. This detection is only possible because we know 
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2000 



Fig. C.2. Power spectra of the expected bias to our lensing 
estimators due to pixelization AC^'^lp.^ (thick solid lines), as 
well as cross-correlation of the pixelization field with the es- 
timators themselves </^^(p^^)^^^ (binned points). The theoretical 

lensing potential power spectrum C^^ is plotted for gradient 
modes as solid black, while the lower dashed black curves give 
(0.1, 0.01, 0.001) X C^^ respectively. The pixelization bias is al- 
ways small compared to both signal and reconstruction noise. 



diM precisely from the satellite pointing. Except at the very high- 
est multipoles, the gradient mode of the deflection lensing field is 
always significantly smaller than the CMB lensing potential. At 
L < 250, which provides most of the weight for the Planck lens- 
ing likelihood, the bias to the lensing potential is significantly 
less 0.1%, while our statistical measurement uncertainty is at the 
4% level. 



information in a multi-million pixel map down to eight eigen- 
modes, it is quite reasonable to expect that the central limit the- 
orem will drive their distribution in this direction. With a set of 
simulated lens reconstructions, it is straightforward to assess the 
Gaussianity of our amplitude estimates. In Fig. D.2 we plot his- 
tograms of the bandpower coefficients from our set of fiducial 
simulations showing that the Gaussian approximation is good 
for the marginal distributions. Further tests of the assumption of 
a Gaussian likelihood are made in Schmittfull et al. (2013) for 
isotropic surveys. There it is shown that parameter constraints 
for a simple two-parameter model of C^j^ scatter across simu- 
lations in a manner consistent with the likelihood width in any 
realization. 

No correlations between the lensing potential and tempera- 
ture power spectra: Lensing may be detected at a significance 
of 1 Oct as a smearing of the acoustic peaks in the temperature 
power spectrum; see Planck Collaboration XVI (2013). A po- 
tential worry when simply combining the temperature likelihood 
with the lensing power spectrum likelihood is that we are double 
counting the same information (Lewis 2005; Perotto et al. 2006; 
Smith et al. 2006) since we do not account for the correlated ef- 
fect of sample variance of the lenses in the estimated C^^ and 
CY spectra. This effect is explored in detail in Schmittfull et al. 
(2013) where it is shown to derive from the connected part of the 
6-point function for the lensed temperature anisotropics. They 
estimate for Planck that the maximum correlation between C^j^ 
and CY for unbinned spectra is around 0.05%, and that ignoring 
this leads to a misestimation of the errors on the lensing ampli- 
tude Al (see Sect. 6.1) in a joint analysis of only 2%. A more 
pressing concern is that the statistical errors on the reconstruc- 
tion power are correlated with those for the temperature power 
spectrum due to sample variance of the CMB anisotropics. This 
effect is also discussed in detail in Schmittfull et al. (2013). The 
correlation derives from the fully disconnected part of the CMB 
6-point function and leads to correlations of at most 0.2% be- 
tween the unbinned spectra. However, these correlations are re- 
moved exactly by the realization-dependent aspect of the N^^"^ 
debiasing procedure discussed in Sect. 2.2. 

To see that the disconnected contribution to cov(C^'^, C^) 
vanishes, we note that with our data-dependent N^^'' subtraction 
the data enters the Ct^ in the form^ 



Appendix D: Likelihood Validation 

In the construction of our lensing likelihood in Sect. 2.3 we have 
used several simplifying approximations to the statistics of the 
reconstructed lensing potential, which we discuss briefly here. 

Uncorrelated binned power spectrum estimates: In Fig. D.l 
we plot Monte-Carlo based estimates of the bandpower corre- 
lation. These plots do not include the efl'ect of beam uncertain- 
ties, errors in the point-source correction or cosmological un- 
certainty in the subtraction of the A/^^^^ bias, which we include 
in our likelihood covariance matrix analytically. Even for adja- 
cent bandpowers, the correlation coefficient is constrained to be 
less than 10% and we ignore such correlations in constructing 
the likelihood. As discussed in Sect. 5, however, we do include 
bandpower correlations from beam uncertainties and the point- 
source and A/^^^^ correction in the likelihood. 

Gaussianity of the binned power spectrum estimates: Our es- 
timator for the lensing potential power spectrum is a sum over 
the trispectrum of the observed CMB. There is no simple, rigor- 
ous argument to expect our bandpower estimates to be Gaussian 
distributed, however keeping in mind that we are collapsing the 



{TiTjTkTi)c{fiTjfkfi 
-\CTlm^CllWj-CTlcil] 

-[ci'TjT,^Cf,TiTi-Cl}c: 



Here, Tf is the inverse-variance filtered temperature, with co- 
variance C^} , and summation over repeated indices is assumed. 
Equation (D.l) arises naturally for optimal trispectrum estima- 
tion of weakly non-Gaussian fields (e.g., Regan et al. 2010). The 
expectation value of the data combination on the right-hand side 
is simply the connected 4-point function, {TiTjTkTi}c, but the 
data-dependent quadratic terms improve the variance of the es- 
timator. The right-hand side of Eq. (D.l) is explicitly symmetric 
under regroupings of pairs of indices so that only the primary 

^ For compactness of presentation, we consider estimates built from 
a single temperature map here. 
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Fig.D.l. Measured correlation coefficients for the bins of 
Table. 1, estimated from 1000 simulated lens reconstructions. 
There is some evidence for correlations along the diagonal at the 
10% level, which is sufficiently small that we have chosen to 
neglect bin-bin correlations in our likelihood. 
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Fig. D.2. Histograms of amplitude coefficients for the bandpow- 
ers of Table. 1, taken from 1000 simulated lens reconstructions. 
Each coefficient has been rescaled by its standard deviation so 
that it can be directly compared to a unit variance Gaussian dis- 
tribution (solid lines). 



coupling of the trispectrum need be included. For lensing, the 
primary coupling factors such that 
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see Hu (2001); Hanson et al. (201 1) for explicit expressions. The 
first set of terms in Eq. (D.l) combine to give 
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M 

where the terms in square brackets are exactly of the form 
of mean-field- subtracted, unnormalized estimates of cpLM and 
(pL-M- The remaining terms are equivalent to the A^^^^ subtraction 
in Eq. (17). Any quadratic estimate of the temperature power 
spectrum will involve a linear combination of terms of the form 
Tp^^Tf\ where T^^^ and Tf^ may be diff'erent noisy estimates of 
the temperature anisotropics but are correlated with themselves 
and Ti through (at least) the common temperature anisotropics. 
If we form the covariance between T^p'^Tf'^ and the data combi- 
nation in Eq. (D.l) that enters the C^j^ estimate, the disconnected 
terms give the contribution to the covariance between the esti- 
mated lensing and temperature power spectra due to sample vari- 
ance of the temperature anisotropics. However, it is straightfor- 
ward to show that by including the additional terms in Eq. (D.l) 
that are quadratic in f the disconnected contribution vanishes. 
We see that, as well as mitigating against errors in the CMB 
and noise covariance matrices (see Sect. 2), and reducing the co- 
variance between C^j^ estimates (Hanson et al. 2011), the data- 
dependent N^^"^ subtraction that we adopt removes much of the 
covariance with the temperature power spectrum. 

The final eff'ect that correlates the lensing and temperature 
power spectrum estimates is the correlation between T and 
due to the late-time ISW eff'ect in the former (Sect. 6.2). This 
produces only local correlations between the spectra and they 
fall rapidly with L being less than 0.5% at L = 40. Their impact 
on the errors on the lensing amplitude Al in a joint estimate from 
the temperature and lensing likelihoods is negligible. 
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